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ANNUAL MEETING OF THE ISRAEL ASTRONAUTICAL SOCIETY 


AN ANALYSIS OF THE LANDING MANEUVER OF WINGED AIRCRAFT 


M. ARENS 
Department of Aeronautical Engineering, Technion-Israel Institute of Technology, Haifa 


ABSTRACT 


The landing maneuver as defined by FAA regulations is analyzed without regard to the 
ground effect. It is shown that if the flare segment is neglected the optimum maneuver 
prior to touch-down consists of a glide segment and one or at most two level deceleration 
segments. If the flare segment is considered, it is shown that for most aircraft the optimum 
landing maneuver consists of a glide segment, a flare flown at maximum Cy, and a level 
deceleration segment. 


1. INTRODUCTION 


Considering the degree of importance currently attached to aircraft landing and 
take-off distance, it seems desirable to analyze the aircraft landing maneuver in 
order to determine the flight path giving the minimum landing distance. Based on 
the requirements of reference 1, the origin of the landing maneuver is defined by 
an altitude of 50 feet, a speed of 1.3 times the stalling speed in the landing confi- 
guration, and an aircraft attitude with respect to the horizontal as determined by 
the steady glide angle at that speed. The landing distance is the horizontal distance 
required to land and come to a complete stop. 

Since the speed at which the aircraft contacts the runway will generally be li- 
mited by safety and/or structural considerations, the landing path may be con- 
veniently broken down into the following two independent segments. 

1. Descent from an altitude f,, velocity V;, and attitude 6; to ground level, 

velocity V2, and zero attitude. 


2. Deceleration from velocity V2 to zero velocity at zero altitude (on the runway). 
Needless to say, minimum landing lengths are obtained when maximum permis- 
sible touch-down speeds are used, due to the large deceleration provided by wheel 
brakes, drag brakes, and thrust reversers. 


Minimizing the first segment requires finding the function @(h) that minimizes 


dh 
| tan 0 


hy 


the integral 
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subject to the equations of motion of the aircraft and the aircraft’s aerodynamic 
and structural limitations. This is a straightforward problem in the calculus of 
variations. The differential equations defining various types of optimium flights 
paths in the vertical plane have been derived by Miele2, however, except for a few 
restricted cases, no general solution can be obtained. However, if the aircraft is 
restricted to small deviations from the glide path, as would generally be the case 
due to safety, aerodynamic, and structural limitations, the optimum flight path 
may be deduced on the basis. of relatively simple considerations. 


The minimizing condition for the second segment is easily found by determining 
the angle of attack at which the maximum braking force is obtained. 


It might be noted that books dealing with aircraft performance generally rest- 
rict their analysis of the descent segment to rectilinear paths>-*. 


2. RECTILINEAR DESCENT 


If the flare distance is neglected and the descent is restricted to rectilinear flight 
paths, the distance to touch down is given by 


OE PE SER UPI E 
pa pri oor ET to eee (1) 
wtf tat ua ee 


hy 1 


If only flight paths of small inclination are considered, so that the lift can be con- 
sidered equal to the airciaft weight, then neglecting ground effects D = D(V), 
since the change in altitude is extremely small. Consideting a constant throttle 
setting (idle) as required by reference 1, the denominator of both integrals is a func- 
tion of V only. Moreover, the second integral depends only on the limits of integ- 
ration, irregardless of the flight path considered. Obviously then, minimizing 1 
requires that the drag minus thrust be maximum during that portion of the flight 
path involving a change in altitude. Thus the optimum descent must be composed 
of a constant velocity glide segment betweer altitudes 4, and zero, plus one or at 
most two constant altitude decelerations, at hy, and/or zero altitude, Figure 1. For 
jet aircraft using flaps the maximum D — 7 will usually occur at V; and the optimum 
descent will consist of a glide to ground level at V; and a deceleration at ground 
level from V; to V2, Figure lc. 


3. FLARE DISTANCE 


For a more accurate estimate of the landing length the necessary change in flight 
attitude should also be considered, even though the “flare” length will generally 
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bea small fraction of the total landing length. Lifting the restriction of rectilinear 
paths to the descent will also change the optimum descent maneuver. 


% 


® 
V =Ny v* vy, 
Fig. 1b Fig. ic 


GOPTIPUI CESCENT (RECTILINEAR) 


Figure 2 
Landing flare 
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Using the notation of Figure 2, the equations of motion of the aircraft are 


pede te 2 

TeWeeD ‘ (2) 
Ww dé 

Se eee 3 

(ay mae (3) 


for aircraft attitudes of small 9. From equations (1) and (2) one obtains 


dV T+ WO — D 
— = (4) 
dl Wig V 


(5) 


where the lift has been expressed in terms of the lift coefficient in equation (5). 


Equation (5) shows that for minimum pull-up lengths changes in attitude should 
be accomplished at maximum permissible lift coefficients and velocities. As will 
be shown later the change in velocity that can be anticipated during the flare is 
small. If the change in velocity is neglected, the shortest flare is one carried out at 
highest permissible C,. 


By introducing the drag coefficient equation (4) can be transformed to 


dV? CpPg 
dl W/S 


V2 = 2g T/W + 2g0 (6) 


The airciaft trajectory, and therefore the horizontal length required for a given 
maneuver, is determined by the simultaneous solution of equations (5) and (6). 
Noting that the term involving © in equation (6) is small and that the velocity 


change during the turn is small, a good approximation can be obtained by sub- 
stituting 


Gi 6, = 25 I TALE dl (7) 
i 
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into equation (6). If attention is restricted to maneuvers flown with C,, equal to 
the maximum permissible C, throughout, equation (6) takes the form 


dV2 cpg 
dl Wis 


8 
V2 = 29(T/W + 91) —2g (5 aad <a)! (8) 
1 


The solution to this equation is 


cDpg 


Bassey! = 
Vie w- |27/s + WIS (20, + ae =) 
Cop Cp 
Cre Cyt. 
ITIS-oW fact ed 2 ee So ares 29 
+ —| |S + W/S (20, + = )| gl ( = ) (9) 


where the subscript | indicates conditions prior to the initiation of the pull-out. 
Substitution of equation (9) into equation (5) gives the differential equation of 
the trajectory. 


dd cypg g 
dl 2W/S = a—bl + ce-? ”) 
where j 
ee w/s (20 2 aaa Sg | 
a =| 271s + W/ ( i+ as ) 


Equation (10) cannot be integrated in closed form. Noting, however, that p << | 
for conventional aircraft, equation (10) can be written as 


ty ety P00 ned ie Bed bi (11) 
dl 2W/S (a+c)—(cp+ d)l 


on 


ce 
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The solution of equation (11) is 


| 2gl Ps 
0,-0 = Aae + ee = taf 4 - 4 (2 -TW- a) | (12) 
: 2(—2 - T/W—8,) he .. ity BS 
Chi 


Computation of some saraple cases indicates that equation (12) gives an excellent 
approximation for the “flare” segment of conventional winged aircraft. In fact, 
in many cases the accuracy obtained by using equation (7) is sufficient. The results 
of a sample calculation are shown in Table [. 


TABLE I 


Sample calculation of horizontal distance, change in velocity and height during flare executed at sea leve! 


WIS TW 0; G Vy & G length AV Ah 
1b/ft2 rad My ft/sec e +s ft. ft/sec ft. 
50 O02? =). 0.1120 ee tieRe 188 2 0.288 186 3 10 


The change in altitude during the pull-up is given by 
Ah = — [tanddl (13) 


Integrating equation (13) one obtains 


a3 2 i 
Pab Seee ogee i ete Pa ATW. 0.) . 
4g(—2- Ne gi ae hy) 


Lt 


inf t = ne = CWS: 0.) * ty ea ee 0) (14) 


ie Cry 


If the time to effect changes in angle of attack are neglected, (J ~ 0), the total des- 
cent distance is given by the summation of the glide, flare, and deceleration distan- 
ces. 

4. OPTIMUM DESCENT 
For an aircraft for which the maximum D— T in level flight occurs at the speed 
V,, the glide and flare should obviously be executed prior to deceleration. However, 
the question arises whether a glide to the height required by the flare, followed by 


Vol. 9C, 1961 M. ARENS 7 


a minimum distance flare and subsequent deceleration represents the optimum 
descent maneuver. Two alternatives may be considered: 
1. Extending the flare portion at the expense of the glide, by flaring at a C, less 
than the maximum permissible value. 
2. Depressing the nose during some portion of the descent so as to accomplish 
a change in altitude at little expense in horizontal distance. 
As for the first alternative, its utility hinges primarily on whether decreasing the 


lift coefficient during the flare will lead to a lower aircraft velocity at the end of the 
flare. 


If for small @, ’s the velocity at the completion of flare is approximated by 


V1= Vee mis (le) 


and the flare plus glide distance by 


Crpg og 2 i i) 
2W/S V2 Wet 
where 
0 
Ah = ———— gh; (17) 
fee aia 
2wis V2 


dV 
It is seen that site 0, and ——<0O except for 


dey, dey, 
‘| 
was 18 
“ ~ ZkWis a) 
it 
Since the lift coefficient required during the landing maneuver exceeds 4kW/S 


for conventional aircraft, this approximate analysis indicates that the flare should 
be flown at the maximum permissible c,. 

The utility of the second alternative, depends on whether the horizontal distance 
required for the dive and subsequent pull-out is smaller than the horizontal distance 
required to glide through the same altitude change. 

Neglecting the velocity changes during the maneuvers involving changes in attitude, 
and assuming each maneuver flown at a constant c, the horizontal length required 
for a maneuver involving a change of attitude AQ is 


AO 
2 ee (19) 
rpg... Gg 


2Wis V2 
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The change in height during the maneuver is then 


A 

kpswmeate ae (20) 
2 

For best results the dive should be executed at minimum c,, the pull-ups at maxi- 

mum c,. For the purpose of this analysis the diving c, will be taken as zero. 


The descent to ground level is thus composed of the following segments, Figure 3: 


8, 


TITIPIPAIT TARE Ae LILLIA TIT 


FIG. 3 
DESCENT INCLUDING DIVE 


1. A diving maneuver at zero c, through the angle A@, followed by a recovery 
at maximum permissible c, through the angle —AQ,. 


2. A glide at the glide angle 0; 
3. A flare at maximum permissible c, to 6 = 0 
The total distance covered neglecting the aircraft’s moment of inertia, will be 


AO? AO Ah 
hanya gr esp en lg ag li CAD 
L ( D -T|W) Po Oe 


2W/S V? 


It is desired to find a stationary condition on /, subject to the requirement that the 
total altitude covered (h;) is constant: 


27 V7 =] 2 
hy = —(Ah),, + {A007 Vin paceman o" pre dig poms 
2 g 71 a J Be al yy 
ZV IS = Ve 2W/S V2 (22) 
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Differentiating eq. (21) we find that a maximum on / is obtained at 
Ad, = 60; (23) 


Therefore, if the nose is depressed through an angle equal to the glide angle the 
landing distance will be a maximum. In fact, any improvement over the conventional 
landing maneuver requires that 


Ad, > 20, (24) 


Since such a maneuver will generally violate safety requirements, and since for 
smaller angle dives the landing distance is increased, it is concluded that the nose 
should not be depressed if minimum landing length is desired. 


GROUND EFFECT 


The analysis of this paper has disregarded the ground effect on the aircraft aero- 
dynamic characteristics. This effect may lead to a significant reduction in the indu- 
ced drag of the aircraft in proximity to the ground. Since the height at which this 
effect becomes significant depends on the wing aspect ratio, wing location, and 
aircraft size, no general conclusions regarding the influence of the ground effect 
on the optimum maneuver can be stated. 


CONCLUSION 


For aircraft for which the ground effect can be neglected, and for which the level 
flight drag minus idle thrust is larger at 1.3 times the stalling speed than at the touch 
down speed, the optimum landing maneuver consists of the following sequence: 
1. Glide at 1.3 times the stalling speed down to the height required for the flare. 
2. Flare at maximum permissible c,. 
3. Level deceleration at ground level to touch-down speed.. 


NOMENCLATURE 
Cp — drag coefficient S — wing area 
Cc, — lift coefficient T — engine thrust 
D — aircraft drag t — time 
g — acceleration of gravity V — aircraft velocity 
h — height of aircraft center of gravity V*— velocity for maximum drag minus thrust 
k — 0¢p/oc? _in level flight 
I — aircraft moment of inertia W— aircraft weight ; 
L — aircraft lift @ — attitude with respect to the horizontal 
i — horizontal distance o — air density 
Subscripts 
1 — origin of landing maneuver f- flare 
2 — touchdown gl— glide 


d — diving maneuver 
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METHODS OF FLIGHT TESTING ABROAD AND IN ISRAEL 


H. Maro 


Flight Test Department, Israel Aircraft Industries 


OBJECT 


The object of this lecture is to summarize briefly methods of flight testing up to 
date, underlining the importance of standardisation and complexity encountered in 
meeting modern engineering requirements and to what extent have we in Israel 
succeeded in keeping up to date. 

In this limited lecture I’d like to bring to the reader’s notice certain of our instru- 
mentation improvements, and the possibilities of applying local know-how using 
already available standard test equipment. 


INTRODUCTION 


On December 17th, 1903, at Kittihawk, USA, Orville Wright flew for 12 seconds 
in a powered aeroplane. Orville and Wilbur Wright utilised just these principles of 
lift stability and control which Leonardo da Vinci had understood 400 years be- 
fore. These men were pioneers in every field—they designed and built their own 
aeroplane, engine and propellers. They even designed their own wind tunnel, they 
were their own stressmen and calculated the loads to which the aeroplane would 
be subjected and perhaps most interesting of all—they were their own test pilots. 

They were the world’s first test team and as such are the finest example for all 
aspiring test teams to learn about and to follow, for they understood their aeroplane 
in every detail. 

There were three main periods of development io flight testing: First when every 
pilot was his own test pilot, then when systematic testing was being introduced 
and finally, when all test flying has been geared to a carefully devised plan. 

It is interesting to note that even with whole batteries of automatic recording 
instruments with telemetering and film monitoring, the test pilot’s humble kneepad 
with its two attached pencils is still as essential as ever, but no longer in the same 
concept as it was understood only a very short time ago when these few scribbled 
down words constituted the feelings, impressions, views, opinions, ideas, halucina- 
ions and even scares of many a test pilot. And even if more than one test pilot 
had the chance to add his own version of the aircraft’s qualities it would again be 
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expressed in a collection of some kind of the above mentioned impressions, which 
when fed back to the design teams brought forth varying results in the final product, 
very much dependent on that individual which up to this time was the min flight 
test instrument, — an instrument of variable reliability, in a science which as the 
number of variable parameters grew was fighting hard to reduce the great number 
of these variables to enable the designer’s equations to have results written out as 
figures and not impressions and feelings. 

It became obvious that the test pilot would have to be eliminated as a reading 
and recording instrument and that not only would he have to be replaced by in- 
strumentation he alone could no longer handle the task of preparing, executing and 
exploiting each flight. 

As in a number of other fields it was the French who soon after World War II 
set up the first Flight Test Personnel School which differed from all existing British 
and U.S. Schools in the fact that not only were test pilots turned out, but the first 
complete test teams trained. Each such team consisted of a Test Pilot, Flight Test 
Engineer, Flight Test Observer and Flight Test Mechanic. 

In practice, to-day it is not the aircraft or the engine which are the subject of a 
specific test program but each quality and quantity of even the smallest modification 
become the subject of project which once submitted for testing by the designers, 
becomes a project of such a complete Flight Test Team. The measurement of 
qualities and quantities which face all those that venture into aircraft design modi- 
fication and even maintenance to-day are common and although in most countries 
they are still called by different names they are more or less derived from the same 
basic exercise of metal moving through the atmosphere. This atmosphere is the 
same for the American Boeing, F-104 and Polaris as it is for the British Comet 
Lightning and Blue Streak and for the French Caravelle, Mirage and Veronique, but 
in each one of these countries, which to-day represent the last word in aviation 
development, was each goal achieved at a different cost to the taxpayer. 

Until recently, and even in many instances to-day, the American and British 
aviation industries have little standardisation in flight testing although every attempt 
is being made to standarise. It is most probably due mainly to France’s limited 
budget that she has succeeded in catching up with the Western World since the war 
but in many a case surpassed their achievements especially in the field of flight test 
instrumentation standardisation. This French success can and is being put down, to 
standardisation and teamwork which are the two main mottos of flight testing 
in France. 

Knowing well the many ways of the many U.S. and British aviation school: 
manufacturers and organisations, each with its own well proven method, it is impor- 
tant to appreciate that in France to-day there is one language in flight testing, one 
catalogue of standard flight test instuments and one manual of flight test methods. 

When the Advisory Group for Aeronautical Research and Development of NATO 
was founded, Theodore Von Karman, its Chairman, wrote a preface which well 
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describes the situation existing, namely, ‘“‘We found there was a great need for a 
flight test manual covering performance, stability and control, and instrumentation 
of aircraft that could be used by design, development or research engineers, test 
pilots and instrumentation personnel of the participating nations in order to expand 
their knowledge, improve their methods, and standardise their techniques. Although 
various members of NATO have their own separate publications covering the 
subjects contained therein, AGARD recognised the need for the compilation, revision, 
and enlargement of this material for the benefit of all NATO nations”. 


STANDARDISATION 


To-day this manual is available to all, not only NATO nations, and with it not only 
the USA and Britain are trying to adopt the standard method as used by France but 
also others like Italy, Germany, Argentina, India and Spain. 


It was not by chence that we in Israel turned for our flight test training to France 
but mainly because we realised that our way if any had to be the cheapest and best 
and most important of all, one which would enable us at all times to turn to someone 
else and compare results. Such comparison of flight test results is only possible 
economically if the two bodies comparing use standard methods and use standard 
recording and measuring equipment. It is at this point that I would like to explain 
what exactly standardisation of test methods and test equipment can do for the 
aviation industy. 


Any potential customer, whether he be in the form of a Buyer, Subcontractor or 
Operator, interested in investigating, modifying or even redesigning his aeroplane, 
is confronted with having to satisfy the requirements spelled out in the language 
of the original airworthiness authority, when he comes to official certification tests. 


It will now save time and expense to a great extent if he is able to carry out hiS 
tests in the same manner, using test equipment knowr to that original authority. 
Knowing which equipment the original designer used his task becomes easier and 
far more re alistic if his results can be compared at every individual stage of testing 
to those of the designer and he can soon spot the point at which he is either erring 
or approaching any red line. If at this stage he is in doubt all he has to do is to advise 
the original designer which pick-ups he used send the standard recorder film in original 
to the designer and he can then await his comments without having to provide him 
with calibrations and reports full of details of the way he measured and trying to 
prove to the original designer that his own non-standard test equipment was not 
at fault in providing unexpected results. 

Let us assume for the purpose of this discussion that some form of standardisation 
has been achieved by most serious concerns and authorities in the world, still Jeaving 
us sufficient elbow room to pick and choose of the various available standard schemes, 
that scheme, which suits the test data acquisition program best economically, time- 
wise and accuracy wise. 
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INSTRUMENTATION 


The basic function of any particular flight test data acquisition program can 
often be satisfied by any of several known schemes, some more widely known and 
used, others less. Each of these schemes, however, has certain short-comings and 
advantages. In order to insure that the basic limitations and capabilities of a particular 
instrumentation system are fully appreciated, it is important that initially the system 
be considered as a whole. 

Instrumentation design under a systems concept requires extensive consideration 
of overall system and cemponent performance, and must be based on a compromise 
of the best features of each subsystem. Rather than adding complexity it is usually a 
requirement to reduce complexity, reduce weight, and increase reliability. System 
approaches to component design often increased system flexibility. This point of 
flexibility will be raised later especially as it concerns our local requirements. 

Through the careful application of these principles, significant advantages can 
be gained in instrumentation systems in terms of increased capability, performance, 
and reduced weight, space and power requirements and improved reliability. 

Another part of flight test data acquisition systems planning is the desirability 
to consider the data reduction and computation process as part of the basic system in 
order to take full advantage of the automatic data reduction aids which have been 
developed. Certain measurement techniques lend themselves more easily to auto- 
matic analysis procedures than others. It is therefore desirable to establish the data 
analysis technique at the same time that an instrumentation system is planned. 


There are eight significant factors which must be considered at the beginning of 
any flight test program before one can begin a detailed plan and design of the instru- 
mentation system let alone the long and usually expensive purchasing and parts 
manufacturing period which cannot be overlooked at the outset. 

Engineering-wise these factors concern both aerodynamic and instrumentation 
facets of the problem and can be summarised as follows: 


1) Number of items to be measured. 
2) The dynamic response required of the measuring system. 


fol 


3) The form of the final record which is desirable. 

4) System complexity which is intolerable. 

5) The degree to which automatic data reduction and computation is both 
feasible and desirable. 

6) The desirability of pre-recording computation. 

7) Required accuracy for the aerodynamic data. 

8) Availability of equipment, facilities and experience. 

In order to indicate the significance of each of the above factors, a short discus- 


sion of each is presented here with an example out of our own short Israeli experience 
where possible. 
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In order to understand the following explanations it should be quite clear to the 
reader that the Flight Test Department receives the demand for flight test results 
in the form of a technical report stating full details of the theoretical work done 
up-to-date indicating which parameters have to be proven in flight, which are required 
for further theoretical work and finally the form and accuracy of the desired infor- 
mation. 


The design and programming of the complete flight test system is now the resposi- 
bility of the Flight Test Team which will now consider the above eight factors to 
determine the best and most economical way of replying to the Engineering Depat- 
ment. 


Number of information channels 


Engineering request demanded the measurement of fuel flows, pressures and tempe- 
ratures, air pressures and temperatures in an existing fuel system in sections where 
modifications were anticipated with the intent of repositioning certain fuel system 
components to make space available for other aircraft equipment. The accuracy 
of the required results was indicated, and obviously, results were required for the 
entire flight domain of the aircraft. Drawings of the system were supplied indicating 
sections where each measurement had to be made. 


ae 


It is easy in any test program to justify any aircraft environmental or meteorolo- 
gical measurement and interestingly enough all these flight variables offer some 
practical interest. There is clearly now, a necessity to determine an order cr priority 
as regards the various parameters it would be desirable to measure. Not only priority 
with respect to sequence or order but certain parameters, will have to be measured at 
the same instant and so the characteristics of the recording systems available will have 
to be examined and due consideration given to the number of parameters required 
at a later stage of the test program. One part of the exercise would be to obtain 
these parameters in the same sections of the fuel system where changes would occur 
in addition to further measurements which have to be made for a new system anyhow. 

The number of variables measured had a lot to do in determining the combination 
weight and size of the instrumentation system. Since the weight and size were some- 
what fixed due to the aircraft structural limitations the only way we could accomodate 
the required number of measurements was by the use of time sharing instrumentation, 
to record on a single channel the large number of temperatures required. On one 
channel we succeeded in getting 16 temperature measurements. This of course 
saved the number of available channels but made the unscrambling more compli- 
cated and brought with it a risk of losing 16 parameters jointly due to the 
malfunctioning of that one channel. 

During this test we used two large A13 Hussenot recorders (see figure 1). giving 
us only twelve traces on each recording instead of a possible sixty which is a theore- 
tical possibility but a film assesors nightmare. 
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Figure 1 
A13 Hussenot Recorder 


Figure 2 


Oscillographs, Quotiometers, Recorder 
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Dynamic response 

The dynamic response required of any measurement system is dictated solely by 
the nature of the experiment and can usually be determined in advance of any testing. 
The dynamic or frequency response of a measurement system is usually established 
by the response characteristics of the transducers and recording devices used in the 
systems. Those transducers and recorders (see Figure 2) which are capable of high 
dynamic responses are also capable of sensing and recording a large variety of 
extraneous noises. Any increase in the signal band width capacity increases the 
noise band width capabilities. 


This is a factor which can cause a large deterioration in the quality of the mea- 
surement record. 


In the particular test of our fuel system the parameters to be measured were 
divided into two groups. The first group consisted of those that vatied slowly; 
temperatures and environmental conditions and the second group varied faster; 
pressures and those prameters which were used to establish the response of the 
flight vehicle itself. As it is impractical to incorporate in a single flight test instru- 
mentation system a variety of recoder reponse characteristics two separate recorders 
(see Figure 3) were used containing recordings of parameters according to their 
qynamic response. 


A2O Recorder 


Figure 3 
A20 and A23 Recorders 
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Figure 4 


Optical oscillograph recording of descent and levelling off in jet A/C. (89 mm film) 
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Record form 


In the matter of the form of the record in which the measurements are to be stored, 
there are two major choices available, visual or electrical. The visual form of record 
can be illustrated by data storage which consists of notes put in a book form by a 
human observer or by lines on a photosensitive film or paper made by an optical 
recording oscillograph (see Figure 4) or by photographs of groups of instruments 
obtained in a photographic observer (see Figure 5). 


Electrical data storage can be accomplished by the creation of a magnetic field 
which is proportional to the quantity being measured, and the storage of this magnetic 
field on a magnetic tape, wire or drum. 


In tests comparing two types of artificial horizons we used a combination of both 
methods. 


A photographictype A20 recorder (see Figure 3) was used to carry all environmental 
parameters as airspeed, altitude, acceleration in Z and X axis, outside air temperature 


Figure 5 


Photographic observer photograph of side of cockpit and position of switches 
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and a time base. A flight test observer was placed behind both instruments which 
were installed at eye level next to each other and he compared the indications against 
a common calibrated scale, noting down all final positions in both form and des- 
cribed the general behaviour of both on a magnetic tape. 

Photography of the instruments was not used as in this sole test requiring a 
cockpit camera installation the method chosen was more economical than installing 
a cinecamera. 

The chief difference between the visual and electrical forms of recording from 
the point of view of the overall system design, lies in the ease with which the data 
can be handled in computation process. 

A visual record is simple and convenient, to determine trends. It must be read 
by a human operator before any further computation can be accomplished (see 
Figure 4). 

An electrical record form is far from flexible and lends itself easily to computation, 
analysis, expansion and mass data consumption which in cases where very large 
aircraft are concerned becomes the solution of the problem. 


System complexity 


An important factor in determining system complexity is the desired scope of the 
instrumentation system. 

If an instrumentation system were planned to solve all of the possible measurement 
problems of the flight tester, there would be produced, no doubt, a system unwork- 
able due to its complexity, size and weight. 

To avoid the creation of such an unwieldy system, the complexity which is permis- 
sible should be clearly established after carefully weighing the advantages and 
disadvantages to be gained by each additional degree of complexity. The factor 
which quite frequently establishes the upper bounds on tolerable complexity is 
often the limitation imposed by the availability, time, funds and personnel 
in a flight test program. 

In one particular test in which the normal position of the aircraft static could 
not be utilised complexity was reduced by applying a principle of comparative 
measurements, An accurate static measurement was made on a standard aircraft 
bya long boom static. Boom readings were recorded for the standard aircraft together 
with its normal static for the entire speed range. When the same boom fitted in the 
same position on the non-standard aircraft was made the comparative recording 
showed the actual difference between the two statics. 


Automatic data reduction 


The choice between manual and automatic data reduction is usually determined by 
the nature of the data and by the amount of computation required, also by the volume 
of data which has to be processed. 
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Up to now most of our data reduction has been manual but with flight test work 
on a prototype in the future, it will become imperative to use automatic data redu- 
ction to keep up with a rapid test program, correcting flaws in between each flight. 
In some cases where the volvme is not large but rather complex, such as in vibration 
measurements, the desired information will have to be extracted by harmonic or 
spectrum analysis techniques. 


Prerecording computation 


We use three methods: mechanical, electromechanical and one that utilizes elec- 
tronic means of performing the computation. 


An airspeed indicator is an example of the first where instead of measuring and 
recording and then calculating from total and static pressures the airspeed, this is 
done by a mechanical computation. An altimeter also falls into this category. An 
illustration of an electro-mechanical type of computation was the use of a trigone- 
metrical resolver to obtain the longitudional and lateral components of a cyclic 
quantity in a helicopter flight test. 


In this type of measurement, cyclicly varying quantities on each of the helicopter 
blades were multiplied by the sine or cosine of the angle between the blade and 
the lateral or longitudional direction. 

The above types of prerecording computation can increase the accuracy, reduce 
the number of channels and reduce materially the amount of data reduction required. 


Prerecording computation has been especially important in Israel when cost 
considerations ruled out the use of automatic data reduction schemes at this 


stage. 


Accuracy 


A major effort in missile and aircraft instrumentation development has been directed 
towards increased accuracy. 

Cost of instrumentation is directly proportional to accuracy to a point where 
from then on accuracy increases in no measurable proportion to cost. Smoothing is 
one of the important methods of increasing accuracy capabilities of a given 
system. 

The importance of a realistic estimate for required accuracy cannot be expressed 
too strongly. The overstating of the required accuracy leads to unnecessary extensive 
labours and needless systern complexity. When the designer requests a degree of 
accuracy without a clear understanding of what this accuracy implies with regard 
to the measurement system, the instrumentation system designer often includes a 
host of nonessential stabilising devices such as power supply regulator, temperature 
compensators etc. which increase bulk and cost of the system. 
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The determination of experimental accuracy has been a source of controversy 
ever since the first measurement was made. The basic reason for this conflict lies in 
our inability to determine true values and the difficulty of one person estimating 
for others the accuracy of his data. Within anyone country, or between countries 
the meaning of accuracy terminology varies, and so here in Israel, where we have 
engineers trained in all four corners of the world the factors affecting accuracy vary 
as much as the different meanings of words like resolution, sensitivity, linearity etc. 


Equipment facilities and experience available 


These factors have the greatest bearing on flight test instrumentation systems. 
Equipment available here ia Israel is mainly French. 

Facilities are of varied origin and our experience is mainly on the equipment 
available in Israel at the moment. 


We have made extensive use of available non standard instrumentation but always 
in such a way as not to affect the final standard recording which is used for compa- 
rison. 

When the choice of a particular instrumentation system is made by assessing each 
of the above factors in accordance with its logical influence on the overall plar, the 
result is a system which is as simple as it can be and yet which satisfies the basic 
aerodynamic requirements. 

The proper balance between instrumentation capabilities and limitations and the 
desirable aerodynamic requirements results in a system which is likely to be reliable 
and compact. 


FLIGHT TECHNIQUE 


Much has been said about instrumentation for the purpose of recording and little 
consideration has been given to the old recorder, the test pilot, and his flight instru- 
ments which have considerable bearing on the accuracy of the entire system. 


Our aircraft have become a complex system requiring accurate handling where 
mechanical controls have been replaced by hydraulic and electric systems. The 
parameters operating them have much narrower margins of operation and require 
to be kept within these narrow ranges of pressures, temperatures and airspeeds 
to keep the entire system functioning properly. 


The test pilot can no longer in addition to flying, keep a log, act as the flight test 
observer etc., but his main job is to fly the aircraft to the limits of human accuracy 
which are with proper training narrower in many cases than those that we are able 
to get from any other type of known flight operator whether remote or automatic. 


Where an average pilot expects to fly within + 100 ft at 40000 feet (see Figure 6a). 
the test pilot is expected, during stabilisation (see Figure 6b), to keep within + 5 
feet, a feat requiring deep concentration and in many cases unusually great effort. 
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Figure 6b 
Recording from test pilot’s stabilisation 
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There is very little point in designing extremely accurate instrumentation when 
the aircraft cannot be flown along predetermined lines to insure the least number 
of variables. 

In the ideal system when flown in perfect stabilisation the only trace which should 
vary should be the one of the parameter under examination (see Figure 6b). 

Apart from a “‘Collimateur” (see Figure 7) the normal instruments of the aircraft 
serve the test pilot in attaining his high degree of flight accuracy. 


LOCAL INNOVATIONS 


An outstanding achievement of local design has been the reduction in the weight 
and the simplification of instrumentation system junction boxes (see Figure 8). 
By using IBM type blocks we have succeeded in reducing the weight and size of 
junction boxes by 90%, a factor in instrumentation of paramount importance, 
especially as the weight of junction boxes has nearly always been an interference 
factor both from the accessibility and adaptability points of view. Using IBM 
blocks, hours of soldering are saved and changes between flights can be made to 
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“Collimateur” sight enabling test pilot accurate flight in X and Y axis. 
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Figure 8 


Junction box using IBM block 


entire circuits within minutes without any changes to microohmic resistances which 
often are the cause of errors in recording. We expect this method to be adopted 
soon universally. 


POSSIBILITIES IN ISRAEL 


We have available in this country both the know-how and instrumentation to enable 
us to undertake a full flight test program of a medium-sized jet transport aircraft. 


The instrumentation, recordings and results attained would be acceptable inter- 
nationally and our results could at any time be compared and cross-checked by 
any authority as all our methods and equipment are standard. 


In addition to basic aircraft design flight test work we can to-day efficiently carry 
out the testing of any modifications and changes to the laid down requirements 
of such authorities as the FAA (U.S.A.), ARB (Great Britain) and others. 
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We can also, by using our equipment and experience, examine existing aircraft 
and equipment in service with airlines and, after carrying out studies of performance 
and operational economy based on accurately measured data, advise in the purcha- 
sing of new equipment. Finally we can test any new equipment considered for pur- 
chase and consult the purchaser how this new equipment will perform in comparison 
to that which he is thinking of abandoning. 


Finally we are available to the Department of Civil Aviation as an arm which 
can be used to determine the performance of any equipment in use for airworthiness 
determination purposes. 


CONCLUSION 


Graphs 1, 2 and 3 show the effect of new flight test methods on the efficiency of 
flight testing. 

The Israel Aircraft Industries, the cradle of Israeli flight testing, have during 
the year 1959-60 test flown @@ types of aircraft and two types of helicopters, reduced 
by 60% the number of flying hours on most types of standard aircraft tested and 
have reduced by 30% the time required to carry out certain tests by introducing 
standard flight test instrumentation. All this with the accompanying reduction 
cost of flight testing where time in the air constitutes the largest single outlay in 
foreign. currency. 

Whatever sucess is claimed to-day by us, or anyone else for that matter, is only 
due to team-work. Team-work of the Flight Test Engineer, Test Pilot, Flight test 
Observer and last but not least the Mechanic, without whose conscientious preparation 
none of this would be possible since the entire flight would be wasted if the aircraft 
would have to return due to unserviceability of one of its systems or, for that 
matter, even any part of the recording system. 
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ABSTRACT 


The addition of nitromethane to various fuels can improve the combustion process. This 
applies, first and foremost, to the decrease of the ignition temperature, and to the reduction 
of the ignition time lag. 

In our experiments, we tested single drops on a hot surface under spheroidal conditions. 
By means of a movie camera, we were able to determine the changes in the drops during 
the evaporation and during the combustion. We also measured the evaporation life of the 
drops at different temperatures. 

In our tests, we used standard fuels which are in common use in engines of various 
types, namely aviation kerosene, diesel fuel, gasoline, and also methy! and ethy] alcohol. 
For the first three fuels mentioned, the nitromethane additions were small (in the range 
of about 2 to 5 per cent) because of low solubility. The solubility of nitromethane in al- 
cohols is unlimited, and these mixtures could therefore be tested for all concentrations. 

We also tested the effect of alcohol additions to nitromethane (under various concen- 
trations) because of the practical importance of this problem, which, at the same time, is 
of theoretical interest. 

Our results showed that small additions of nitromethane produce a great improvement 
on the ignition of fuels. A nitromethane addition of from 0.5 to 3 per cent reduced the 
ignition point and the ignition time lag appreciably under the prevailing experimental 
conditions. At high nitromethane concentrations in alcohol further additions of nitro- 
methane can, in a certain range, have no further effect on the ignition, and can even produce 
a negative effect. For alcohol-nitromethane mixtures, the ignition point can be lower than 
for either constituent alone. 

Measurements of the flashpoints (Martens-Pensky method) showed that the addition 
of nitromethane to fuels (at low concentrations) causes reductions that are not appreciable, 
although the flashpoints of the mixtures can be lower than those of their constituents 
alone. 

The effect of small nitromethane additions to fuels on the evaporation time is not great, 
but such additions produce changes in the “kinematics of evaporation”’ in the sense that 
they change the evaporation rate of the given drip at certain stages of its evaporation 
life. This phenomenon seems to be connected with the effect on the ignition. 

The results of the tests can be important from the practical point of view for different 


types of combustion engines using fuel injection. 


I. INTRODUCTION 


Various investigators have, for some time, tested the physical and chemical proper- 
ties of nitromethane!. As is well known, it is difficult to use nitromethane as a sole 
engine fuel, first and foremost because of the great danger of explosion. However, 
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nitromethane can be utilized as one component dissolved in various other liquid 
fuels for either jet or piston engines. One can also use small quantities of nitromethane 
as additive to various liquid fuels in order to improve those properties of the fuels 
that are most important in actual engine operation. Although experiments along 
this direction have been carried out and a number of patents recorded 1.2.3, the 
problem has not yet been sufficiently studied, and, practically speaking, has not 
as yet emerged from the walls of the laboratory. 

The addition of nitromethane to various liquid fuels can improve the combustion 
process4. By means of such additions, it is possible to reduce the ignition point, and 
also the to reduce the ignition time lag and the combustion time. The ignition pro- 
blem is closely connected with the problem of drop evaporation. This is, of course, 
generally true, and not only for nitromethane additions. 

The solubility of nitromethane in a number of fuels is small. In our measurements 
we obtained the following solubilities at 20°C, in percent (by weight): 


In diesel fuel: 2.1 percent 
In kerosene: 3.2 percent 
In gasoline: 5.7 percent 


The solubility of nitromethane in ethyl or methyl alcohol is unlimited. 


Our measurements showed that the reduction of the fuel flashpoints that are 
caused by the nitromethane additions are not appreciable, so that, practically spea- 
king, there should be no problem in using such additions (Figure 7). 


II, EXPERIMENTAL SET-UP AND PROCEDURE 


In our investigation, we tested single drops on a hot surface. The following variables 
were measured: 
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where t is the evaporation time of a single drop in seconds, f, is the temperature of the 
hot surface in °C, t,, is the minimum hot-surface temperature necessary for ignition 
o is the precent (by weight) concentration of nitromethane in the fuel tested, and 
At,, is the ignition time lag in seconds. 


The tests were carried out as follows: Using a special dropper (Figure 1), drops 
of variable size vcan be obtained. The upper end of the dropper is connected with 
a long rubber tube passing between two hard rubber rollers. By rotating roller 5 
(Figure 1), the section of the tube between the dropper and the rollers is lengthened 
or shortened, thus increasing or reducing the minimum pressure required for relea- 
sing a drop. This results in drops of uniform shape and size, irrespective of the 
height of the liquid column in the dropper. 
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Figure 1 


: Drop-release instrument 
i—Dropper; 2—Stand; 3—Turning mechanism; 4—Handle; 5—Drop-squeezer; 6—Shielding door; 
7—Plate. 


The set-up is such that the experiment can be carried out from behind a shield. 
Thus, the dropper is charged and manipulated from behind the shield by means of 
a suitable mechanism, so as to prevent heating of the dropper and of the liquid in it 
(both for reasons of safety and of accuracy), and to return it to its correct position 
over the center of the plate for each drop. The shield is provided with a steel door 
having a plexiglass observation window. The set-up is shown in Figure 2. 


Experimental set-up 
1—Heating plate; 2—Base; 3—Dropper; 4—Pipe rollers; 5—Rubber tube; 6—Handle; 7—Stand; 
8—Thermocouple; 9—Galvanometer; 10—Thermometer; 11—Transformer; 12—Variac; 13— 
Cooling water; 14—Cooling pipe; 15—Movie camera. 


The evaporation time t was measured with a stop watch, and also by means of 
photographs taken with a movie camera. The plate is made of stainless steel. The 
contact surface must be very smooth, and requires frequent polishing in the course 
of the experiment. The surface temperature ¢, was measured by means of a chromel- 
aiumel thermocouple attached to the plate with a screw. 

The above-described method for studying evaporation has already been used by 
other investigators®. 

The movie film (taken at 64 frames per second) shows the changes taking place 
in the drops while on the hot surface, and also enables us to measure the ignition 


time lag At,,. 
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The evaporation curve t = f; (¢,) can be divided into three sections (Figure 3). In 
the first section, where the temperatures are low, an increase in surface temperature 
t, Sharply reduces the evaporation time t, up to a certain minimum value. After 
this minimum, the second section begins, in which an increase in f, sharply 


200 50 800 y 400 800 
Surface temp,t,[‘c] Surface temp, t,- [°c] 
50 a 


200 500 800 0 400 800 
Surface temp. t-[c] Surface temp. t- [ec] 


AY 
Figure 3 


Evaporation time t (seconds) of drop on hot surface as function of surface temperature ts (°C) 
(Drop diameter approx. 3 mm.) 
Curves A: a—kerosene; b—kerosene with 1.5 per cent nitromethane. 
Curves B: a—diesel fuel; b—diesel fuel with 1 per cent nitromethane. 
Curves C: a—gasoline; b—gasoline with 3.0 per cent nitromethane. 
Curves D: a—methyl alcohol; b—methy] alcohol with 10 per cent nitromethane 
The ignition points are marked on all the curves. 
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increases t up to a certain maximum value. After this maximum, the third section 
begins, in which an increase in f, again reduces t. 

The first section of the curve represents evaporation under non-spheroidal condi- 
tions. (The drop loses its spheroidal shape immediately upon coming into contact 
with the hot surface, and is transformed into a thin film covering the whole hot 
surface.) In the third section of the curve, the drop retains its spheroidal shape until 
it is completely evaporated, which represents evaporation under spheroidal condi- 
tions4, The second section represents a transition between the conditions prevailing 
in the first and the third sections. 

In our experimental method, the values of the ignition temperature represert the 
temperatures ¢;, of the hot surfaces upon which the drop is placed. The experiment 
was carried out under ordinary ambient temperature and pressure conditions. The 
measurements of the ignition time lag At,, were carried out by means of the same 
method and under the same conditions. 


In additional experiments, the flash points of methyl alcohol, of ethyl alcohol, and 
of kerosene, with added nitromethane were measured. These measurements were 
carried out by means of the Martens-Pensky apparatus, with initial cooling. 


III. RESULTS OF THE MEASUREMENTS 


Figure 3 shows the results of the evaporation-time measurements (single drops on 
a hot surface) for the following fuels: 


Curves A: Aviation kerosene 

Curves B: Diesel fuel 

Curves C: Commerical automobile gasoline 
Curves D: Methyl alcohol (99 percent) 


The curves “a” and “‘b” represent the various fuels without and with nitromethane 
additions respectively. The following percent (by weight) nitromethane additions 
were used: 


1.5 percent in kerosene 
1 percent in diesel fuel 
3 percent in gasoline 
10 percent in methyl alcohol 


On each of the above curves, the ignition temperature (in °C) is also marked. The 
diameter of the drops was approximately 3 mm. 

The curves of Figure 3 show that the effect of the nitromethane additions on the 
evaporation time is not great. The additions to kerosene (A) and to diesel fuel (B) 
increase the evaporation time t in the spherical range, while they decrease 1 for 
gasoline (C) and methyl alcohol (D). In all cases, the nitromethane addition appre- 
ciably reduces the ignition temperature (Figure 3). 
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Figure 4 presents results of our measurements of the ignition temperature 14; as 
a function of the nitromethane concentration o in industrial ethyl alcohol - curve 
“a” and in methyl alcohol (99% grade) -—curve “‘b”. The curves show that an 
increase of the nitromethane concentration in these alcohols reduces the ignition 
temperature ¢,, only within a certain range, up to o © 25 percent (by weight) in 
our experiments. For o > 25 percent, further increases in o produce either no 
appreciable changes at all, or a slight increase in ¢,,. It is characteristic that, for o 
close to 100 percent, a further increase in o increases ¢,,. This means that, if we 
add methyl or ethyl alcohol to nitromethane, we can in this way also reduce the 
ignition temperature t,,. 


We thus see from Figure 4 that it is possible to obtain 


(tig) 


and also (tip) 


(tig) ot 


(tig) 3 


Fi+F2 


Fi+F2 


where (tig) 4 is the ignition temperature for the first fuel, 
(tig) 3 is the ignition temperature for the second fuel, 
and On is the ignition temperature for the mixture of the two fuels. 
In other words, the ignition temperature of fuel mixtures can be lower than the 


ignition temperature of each component separately. 


500 


@-PER-CENY (BY WEIGHT) 


Figure 4 


Ignition temperature fig (°C) as function of nitromethane concentration in per cent (by weight) 
in ethyl alcohol—curve “‘a’”, and in methyl alcohol—curve “‘b’’. 
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@~PER-CENT (EY WEIGHT) 


Figure 5 


Ignition temperature fig (°C) as function of ethyl alcohol concentration ¢—curve ‘‘a’’, and methyl 
alcohol concentration @—curve “‘b”, in per cent (by weight) in nitromethane. 


©-PER-CENT; (BY WEIGHT) 


Figure 6 
Ignition temperature fig (°C) as function of nitromethane concentration @ in per cent (by weight 
in diesel fuel—curve ‘“‘a’’, kerosene—curve ‘“‘b”, and gasoline—curve ‘‘c’’. 
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o-PER CENT (BY WEIGHT) 


Figure 7 


Flash point (°C) of ethyl alecohol—curve “‘a”’, methy! alcohol—curve “‘b’”’, and kerosene (aviation)— 
curve ‘“‘c’”? as function of nitromethane concentration o in per cent (by weight). 


It might be worth nothing that Barrére and MoutetS have pointed out the possibi- 
lity of obtaining an ignition time lag of a fuel mixture which is smaller than the 
corresponding Jag for each individual component. 


Figure 5 presents results of our measurements for additions of ethyl alcohol 
(industrial)—curve ‘‘a’”’, and additions of methyl alcohol (99 percent grade)—curve 
“b”, to nitromethane. The curves show that additions of these alcohols reduce the 
ignition temperature ¢,, of nitromethane. This phenomenon can obviously be of 
practical importance. 


Figure 6 shows results of our measurements for additions of nitromethane (in 
small quantities) to diesel fuel— curve “a”, to aviation kerosene — curve ‘“‘b’’, 
and to commercial gasoline — curve “‘c”. The curves show that the nitromethane 
addition has an appreciable effect on the ignition temperatures ti, of these fuels. 


Table I shows the effect of nitromethane additions on the ignition time lag Ar,, 
of kerosene, as found in our measurements carried out by means of a movie camera 
photographing 64 frames per second. The temperature of the hot surface had to be 
adjusted according to the respective ignition temperature t;, for the given mixture. 
Table I shows results for two different values of hot-surface temperature. 
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TABLE I 
Ignition time lag ATig of single drops on a hot surafce of temperature ts 


Nitromethane concentration o per cent 


(by weight) in aviation kerosene 0 0.5 1.0 3.0 
Temperature of the hot surface fs (°C) 675 675 552 552 
Ignition time lag ATig (seconds) 0.36 0.09 0.30 0.17 


Figure 7 shows the effect of nitromethane additions on the flashpoint of ethyl 
alcobol — curve ‘‘a’’, methyl alcohol — curve “‘b’’, and kerosene— curve “‘c”. The 
additions produce a reduction of the flashpoint. 

By means of movie-camera photographs of single drops on the hot surface, the 
changes taking place in the drop during evaporation, ignition, and combustion can 
be studied. For the above-mentioned fuels used in this investigation, and also for 
some other fuels, we have established the fact that the drop at certain moments 
experiences an appreciable increase in diameter, while, at other moments, the drop 
diameter suddenly shrinks. This change in the drop diameter is liable to appear at the 
beginning, and also several times afterwards in the course of the evaporation and 
combustion process. For diesel fuel, these changes in the drop diameter are es- 
pecially pronounced. When nitromethane additions are used, this phenomenon 
appears to an even greater extent. 

These changes undoubtedly influence the evaporation end ignition processes 
and the combustion time, but it is difficult to determine any systematic correla- 
tion. It seems that the phenomenon might be of importénce in trying to explain 
the effects of additions to fuels on their evaporation and ignition, and also the 
connection between the evaporation and ignition of drops of various fuels. 


IV. CONCLUSIONS 


The addition of nitomethane to various fuels (diesel fuel, kerosene, gasoline, methyl 
and ethyl alcohol) influences their evaporation and ignition processes. 
1. Effect of nitromethane additions on the evaporation of drops. 

a. For some fuels (diesel fuel and kerosene), the nitromethane additions cause 
an increase in the evaporation time, while, for other fuels (gasoline and methyl 
alcohol), the evaporation time is decreased. The effect, however, is not appreciable, 
and, in our tests, amounted to approximately 10 percent. 

b. During the evaporation process, vapours (and also gases resulting from pos- 
sible decomposition) can collect inside the drop. The drop expands up to a certain 
size, whereupon it can burst, with the consequent release of the vapours or gases. 
This produces a new drop of smaller diameter taking the place of the previous drop. 
The new drop again expands, and this sequence can be repeated several times until 
the end of the evaporation process. This phenomenon is very pronounced for diesel 
fuel, and becomes even stronger with nitromethane additions. 
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2. Effect of nitromethane additions on the ignition temperature ¢,, and on the 
ignition time lag Ar,,. 

a. Small concentrations of nitromethane cause a great decrease of 1,,, especially 
for diesel fuel, kerosene and gasoline. Thus, 0.5 percent (by weight) of nitromethane 
in kerosene or in diesel fuel reduces ¢,, by about 100°C. In order to obtain such a 
reduction in alcohols, approximately 10 percent of nitromethane is required for 
ethyl alcohol, and approximately 20 percent in methyl alcohol. 

b. The ignition temperature ¢,, for a mixture of two fuels can be lower than that 
of each fuel alone. Thus, the ¢,, of ethyl alcohol and that of nitromethane were, in 
our tests, found to be 640°C and 535°C respectively, whereas the minimum 1,,for a 
mixture of these two fuels was 500°C. 

c. The ignition time lag At,, of the various fuels tested is reduced by the addition 

of nitromethane. Thus, a 0.5 percent addition of nitromethane reduces At;, of 
kerosene by a factor of 4. 
3. The addition of nitromethane to various fuels reduces their flashpoint, but 
by a relatively small amount. Thus, a 0.5 percent addition of nitromethane to kero- 
sene reduces its flashpoint by 7°C. It is characteristic that the flashpoint of a mixture 
of two fuels also can be lower than that of each fuel alone. This has been verified 
for mixtures of nitromethane with methyl and ethyl alcohol and also with kerosene. 
4. The addition of methy! or ethyl alcohol to nitromethane of up to about 20 per- 
cent reduces its ignition temperature ¢,, by as much as 30°C. 
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STRUCTURE, VARIATIONS AND MEASUREMENTS OF THE EARTH’S 
IONOSPHERE AND EXOSPHERE 


COLMAN ALTMAN 
Department of Physics, Technion — Israel Institute of Technology, Haifa 


ABSTRACT 


The extensive use of high altitude rockets during, and subsequent to, the IGY period has 
yielded much useful information about the structure of the E and F layers of the ionosphere. 
Measurements in earth satellites and space probes have revealed the inner and outer rad- 
iation belt, and also ring-currents encircling the earth at distances of many earth radii. 
These measurements have supplemented and enlarged the existing body of knowledge of 
the inner ionosphere (D, £ and F regions) acquired primarily by ionosonde techniques 
and geomagnetic measurements, and also of the outer ionosphere (or exosphere), in which 
atmospheric ‘“‘whistlers’’ played a pioneer role. Our present state of knowledge is surveyed 
and work being conducted in this field by the Ionospheric Physics group in the Technion 
is briefly described. 


INTRODUCTION 


The movement of the earth and its dipole magnetic field through a very hot hydrogen 
gas — an extension of the sun’s coronal atmosphere — and the consequent conduc- 
tion of heat into the earth’s atmosphere plays an important role in the thermal! 
economy of the high atmosphere. The ‘solar wind’ consisting of neutral couds of 
ionized hydrogen is effective in sweeping away the extension of the earth’s magnetic 
field at a distance of some 5-8 earth radii on the sunlit side of the earth and creating 
a magnetic discontinuity in part of the auroral zones into which high speed protons 
are funnelled directly from the solar wind. In addition, high speed protons and 
electrons which leak through the magnetic barrier are trapped in the geomagnetic 
field and give rise to the outer Van Allen radiation belt. This belt, together with the 
direct contribution of solar protons in the near polar regions, heat the auroral 
atmosphere above 100 km height and generate wind circulations which transport 
beat to the equatorial and temperate zone ionospheres, and heat the thermosphere 
of the earth to over 1000°K. 

The ultraviolet radiation of the sun and soft X-rays ionize the D-, E- and F- 
layers of the ionosphere, and the interplay of production, collision and recombination 
rates of ion pairs generate the variation of electron densities in the ionosphere which 
determine the conditions of electromagnetic wave propagation. Diurnal tides in the 
ionosphere generated by solar heating institute world-wide current systems in the 
E-layer which are manifested in the diurnal geomagnetic variation. 
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At the base of the thermosphere water-vapour and methane are decomposed by 
solar ultraviolet radiation and by upward diffusion the neutral hydrogen atoms 
from a “telluric hydrogen corona”. Thermal protons are generated by charge- 
exchange reactions snd their relative abundance gradually increases until they 
become the main component at a height of one earth radius. 

The Van Allen radiations, formed by trapped solar emissions, although important 
in their effects, form a negligible proportion of the total ion content at these heights, 


Each component — the neutral hydrogen atoms, the thermal protons and the 
energetic Van Allen radiations —is manifested by its distinctive effects. Some effects 
and proposed measurements are discussed. 


THE INTERPLANETARY MEDIUM 


The medium in which the earth’s orbit issituated is determined by conditions pre- 
vailing in the sun, constituting a distant extension of the solar corona. At a height 
of 50,000 km above the solar surface, the corona acquires a kinetic temperature of 
one to two million degrees, consisting mainly of a completely ionized hydrogen 
plasma. The thermal conductivity of a plasma is proportional roughly to 75/2 and 
because of the high conductivity attained, the corona is comparatively isothermal 
over large radial distances. At the earth’s orbit Chapman!.2 finds that this plasma 
has a temperature between 50,000°K and 200,000°K (corresponding to 6-25 ev 
protons), using a model of hydrostatic equilibrium in which T = 106°K for the 
corona in the sun’s vicinity. Using higher values of coronal temperature — 2 x 106°K 
— which certainly exists in active regions of the corona, Parker3 has found that 
the hydrostatic equation is not valid, yielding finite values of density at an infinite 
distance, and concludes that the corona is continually expanding outwards, as a 
“solar wind’. Using an equation of hydrodynamic equilibrium Parker finds solar 
wind velocities of the order of 500-1500 km/sec at the earth’s orbit (up to 10 kev 
proton energies). 

Experimental observations support these conclusions. Biermann45 has found 
that extreme accelerations of comet tails can only be explained by corpuscular 
streams from the sun of densities between 50 and 104 particles/cm3 and velocities 
of some 103 km/sec at the earth’s orbit. Measurements on the zodiacal light® have 
yielded electron densities at the earth’s orbit of 50-104 electrons/cm3. The figures 
vary from low values at sunspot minimum to high values in periods of intense solar 
activity. 

Finally we should mention the close correlation between solar flares and geomag- 
netic activity on the earth 24-40 hours later, which clearly indicate that the ionized 
clouds emitted by the sun travel with radial velocities between 1000-1500 km/sec. 

To summarize, the picture which emerges is roughly as follows. The earth’s orbit 
is immersed in @ hot solar ionized gas of temperature between 50,000°K to 200,000°K 
with a density of several hundred protons/cm3 through which a solar wind of some- 
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what lower density frequently blows, tending to carry with it the stationary inter- 
planetary gas. At times of solar activity, the solar wind takes the form of dense, 
high velocity ionized couds. Because of its relatively high kinetic energy density 
the ionized solar wind will tend to sweep away all transverse magnetic fields in its 
path, leaving behind it radial magnetic field lines, which slowly spiral outwards 
because of the sun’s rotation. 


THE INTERACTION BETWEEN THE SOLAR WIND AND THE EARTH’S MAGNETIC FIELD 


Equating the magnetic energy density B2/8x of the earth’s dipole field with the 
kinetic energy density 4} n m v2, we find equality at radial distances of 5-8 earth 
radii on the sunlit side. We may thus expect the earth’s dipole field to be swept 
away at larger radial distances, whereas at small distances the solar wind wil] be 
excluded. We thus define an outer limit to the outer radiation belt. If penetration 
by high speed protons occurs, these particles will be trapped and constrained to 
spiral back and forth between magnetic mirror points along lines of magnetic field. 
The outer limit of the geomagnetic field will correspond to field lines which meet 
the earth’s surface at geomagnetic latitudes of 65°-70°, i.e. in the auroral zones. 
Furthermore, at somewhat higher latitudes (towards the polar regions) clouds of 
solar protons will frequently pour into the earth’s atmosphere on the sunlit side of 
the earth. Such clouds probably generate by travelling-wave-tube mechanisms? 
or cyclotron radiation’ the audio frequency electromagnetic emissions such as 
““dawn-chorus’’, hiss, etc., which in high latitudes are regularly received in the early 
morning and in times of geomagnetic activity7.9»10, 

On the night-side of the earth the magnetic field lines are “‘shielded” from the 
solar wind, and hence will extend outwards considerably, some investigators belie- 
ving up to one hundred earth radii. 

A further postulated effect of the solar wind interacticn with the geomagnetic field 
is that electrons and protons being oppositely deflected will set up electrostatic 
forces which will drive “ring currents” to girdle the earth!2. Magnetic measurements 
in the first Soviet moon rocket!3 revealed such a westward flowing ring current 
at a height of 22,000 km from the center of the earth (inside the outer Van Allen 
belt), and later measurements 14 in Explorer VI revealed an additional (or perhaps 
displaced) ring current system of some five million amperes girdling the earth at 
a radial distance of 40,000 km. 


THE RADIATION BELTS AND THE EARTH’S EXOSPHERE 
The outer radiation belt particles are engendered by penetration of electrons and 
protons having an energy peak probably near 10 kev. Once captured, they spiral 
back and forth along lines of magnetic force, their pitch angle 6 (angle between 
guiding line of force and direction of velocity vector) being related to the magnetic 
intensity, B, by the relation 
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B/sin2 8 = const, and the mirror point, M, being defined by 0,, = 7/2, and 
B,, = B/sin20. The steeper(smaller) the initial angle, 0, the closer to the earth’s 
surface will the mirror point be, and the smaller the survival time of such a particle. 
Hence a non-isotropic distribution is soon set up in which large pitch angles are 
favored. Energetic particles are lost to the medium following charge-exchange col- 
lisions with neutral hydrogen atoms in the belt. The collision rate for such collisions 
has been shown to be constant for energies between 100 ev and 20 kev above which 
the collision rate falls rapidly!5. We may thus expect the outer Van Allen belt radia- 
tion to contain protons with energies well above 20 kev, and according to Stuart! 
with energies mainly above 50 kev. The electron component has a much longer life- 
time and is the main radiation constituent of the outer belt, concentrating in the 
energy range 20-100 kev!7. Nevertheless their absolute abundance is very low, 
roughly 10-! electrons/cm3 with energies greater than 20 kev18. 


The generation of the outer belt radiation in the solar wind has been clearly demor- 
strated. Measurements of the drag on Sputnik HI!9 (which passed regularly through 
the outer belt in the auroral zone) showed a 27-day cycle and also two large increases 
in drag, each time some 36 hours after major flares in the sun, and in step with 
geomagnetic storms. Vanguard I, on the other hand, which passed through the 
inner radiation belt showed little such effect. This lent support to the theory that 
the inner belt originates primarily in the decay of neutrons in the cosmic ray albedo, 
which would account too for the relative hardness of its radiation. More convin. ing 
were the results of the Pioneer IV flight!8, launched six days after a strong magnetic 
storm, which revealed radiation intensities about 100 times higher than previous 
measurements by Pioneer [If and by Mechta, the Soviet moon rocket, during mag- 


netically quiet periods, indicating clearly injection of corpuscular radiation by 
ionized solar clouds. 


THE THERMOSPHERE AND HETEROSPHERE OF THE EARTH 


The temperature of the atmosphere reaches a minimum (about 150°K) at a height 
of some 85 km in the so-called mesopause, and above this the temperature gradient 
is positive—definining the earth’s thermosphere. Nicolet and Chapman consider 
that the thermosphere continues (i.e. that its temperature gradient is positive) until 
it finally merges with the hot solar plasma beyond the radiation belts2°, Rocket 
measurements however, based on pressure-altitude variations21!,22, indicate that the 
atmosphere becomes roughly isothermal close to 200 km with temperatures varying 
from 1000°K-2000°K depending on local time and season. At an altitude of about 
550 km the mean free path in a horizontal direction of the neutral atmospheric 
components equals the scale height, so that neutral particles travelling vertically 
with sufficient energy are able to leave the earth’s field without further collisions. 
This leve: is the base of the neutral particle exosphere, in which temperature is no 
longer well defined, end in which usual gas-kinetic relationships no longer hold23. 
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Above the mesopause solar ultraviolet radiation is strongly sbsorbed, dissociation 
of Og and to a lesser extent Np resulting. Because of the low recombination rates, 
which decrease with altitude, we find that the concentration of atomic oxygen in- 
creases and it becomes the main atmospheric constituent at a height of several hund- 
red kilometres, although at 250 km molecular nitrogen, No, is probably still more 
abundant?4, 

Because at these altitudes vertical diffusion is much more important than mixing 
due to wind circulations, the density of each constituent tends to fall off with altitude 
according to its own scale height H = RT/Mg (i.e. diffusive equilibrium is established) 
and the composition of the atmosphere changes constantly with height, lighter ele- 
ments gradually predominating. This region is thus called the earth’s heterosphere. 
The decomposition of water vapour and methane just below the mesopause2" liberates 
neutral atomic hydrogen which, because of its large scale height (even though scale- 
height is not strictly defined in an exosphere), replaces oxygen as the main neutral 
atmospheric constituent above 1200 km25, and forms a “‘teJluric hydrogen corona” 
at greater heights. The presence of this vast hydroger cloud or corona encompassing 
the earth (with particle densities of 105 particles/cm3 at 200 km and 3 x 104 part- 
icles/em3 at 1200 km25) has been demorstrated by rocket measurements of the 
night sky26 which was found to be aglow with diffuse Lyman-a« emission (rescnantly 
absorbed and re-radiated from the Lyman-a« flux of hydrogen in space), and also 
by the shape of the narrow absorption strip in the centre of the solar Lyman-« 
emission line27. 

The atmospheric densities at heights ranging from 100 km to 700 km have been 
measured by calculation of drag on sateilites!9, and directly by rocket-borne ioni- 
zation gauges2!. These measurements were interpreted (for instance by Jastrow28) 
as indicating the simultaneous occurrence of temperatures of 2000°K in the auroral 
zones and 1000°K in the temperate-equatorial zones at heights of 200 km, corres- 
ponding pressures in the auroral zones being five times higher than in the temperate. 
Johnson later demonstrated that such pressure and temperature differences were 
dynamically inconceivable at these heights?9. For even under extreme anticyclonic 
conditions (which apparently do not exist at this altitude) the constant pressure 
surface would be 12 km higber over the poles than over the equator. Additional 
pressure differences must introduce a wind circulation from pole to equator. Now 
heating of the auroral atmosphere by energetic radiations in the outer Van Allen 
belt and solar corpuscular streams spiralling in along field lines outside of the outer 
belt undoubtedly does occur. Heating would occur mainly at 100 km altitude tailing 
off towards 200 km as evidenced by the height distribution of auroral displays. 
We may thus assume (with Johnson?) a circulation in which high wind velocities 
would be set up at 200 km and above (dropping with a steady gradient to zero at 
100 km), which would transport heat from the polar to the equatorial atmosphere 
and here downwards from the observed temperature of 1500°K at 200 km to 200°K 
in the mesopause (~ 100 km). Assuming a wind velocity of 150 km/sec at 200 km 
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and above we may show that sufficient heat would be transported to supply all the 
heat requirements of the temperate and equatorial ionospheres. Further a circulation 
of this type would be maintained by a constant pressure surface only 8 km higher 
at the pole than at the equator at 200 km. Adding this to the maximum 12 km 
height difference due to conceivable limiting anticyclonic circulations, and bearing 
in mind that the scale height at this level is some 50 km, we find that a 50% pressure 
difference (which could be engendered by a 100°K temperature difference only 
between poles and equator) would set up a wind circulation which would cover 
the thermal budget of the entire ionosphere. It would appear now that the extreme 
pressure and temperature differences reported stem from the fact that the auroral 
measurements were performed during a more active phase of the sunspot cycle than 
the corresponding temperate zone measurements, a factor which presumably had a 
considerable influence. 


At higher altitudes (at 400 km say) the auroral pressures might indeed exceed the 
temperate values by a factor of five, since the viscous foices in a given wind velocity 
gradient are independent of pressure, and so become a dominant factor at higher 
altitudes. 


To conclude this section it may be remarked that measuremenst of wind velocities 
above 90 km by radio observations of drifting meteor trails3° have shown that above 
105 km winds blow from pole to equator all year round, summer velocities reaching 
40 metres/sec at 105 km, and the velocity gradient is always positive above 90 km. 
This seems to be clear support for Johnson’s model. 


THE IONOSPHERE OF THE EARTH 


The ionosphere is formally considered to start where the ionization density is suf- 
ficiently high to affect radio propagation, and this is just below the mesopause 
(~ 80 km). For many years the ionosphere was investigated by ionosonde techniques 
only, which yielded fairly accurate height versus electron-density profiles up to the 
level of maximum ionization. These profiles tended to show distinct layers of ioni- 
zation termed the D, E and F layers according to the height at which they occurred. 
However, recent rocket measurements have shown that frequently “ledges” (where 
the density gradient is positive but small) tend to form between layers rather than 
“valleys”. In view of this a formal definition has been adopted in which the D-layer 
is that part of the ionosphere below 90 km, the E-layer is situated between 90 and 
160 km, and above 160 km the F layer is located. 


The F-layer is formed primarily by the ionization of atomic oxygen by the entire 
spectral range of solar radiation from 200-850A,32 so that ionized atomic oxygen 
is already the main ionic component of the atmosphere at 200 km, as verified by 
rocket measurements33. The E-layer is formed mainly by soft X-rays in the 10- 
100 A range?? and the D-layer by harder X-rays (A < 10A) and, proceeding to 
lower altitudes by Lyman-« and by cosmic radiation34. It should be remembered, 
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however, that these ionospheric layers are predominantly neutral regions. The ion 
to neutral component ratio at the E and F layer peaks are about 10-7 and 10-4 
respectively. 


Above the F-layer peak recombination gradually becomes unimportant and a 
diffusive equilibrium is set up for the ionic component too. Now hydrogen ions 
are formed by charge-exchange between hydrogen atoms and oxygen ions, the ioni- 
zation potential being almost identical for the two35. Charge exchange collisions 
will be important only up to the base of the exosphere, and since the reaction can 
proceed equally well in both directions the neutral and ionized oxygen-hydrogen 
ratios here will be equal, and relative ionic scale heights will determine the subsequent 
ionic distribution higher up. We thus find that at 1600 km protons become the main 
ionic component25, 


Now the scale height of a plasma is determined by a temperature equal to the 
sum of the ionic and electronic temperatures, so that the ionized hydrogen scale 
height is about twice that of the neutral hydrogen. (Note that because of the large 
cross-section for ion-ion collisions, the ion-ion mean free path is nowhere larger 
than the scale height, and the ionized hydrogen component is nowhere in an ionic 
exosphere, but instead governed by hydrostatic equilibrium). A second factor which 
increases the scale height of the ionic component is the magnetic field. Because of 
the high conductivity of the ionosphere magnetic field lines are “anchored”’ in the 
lower ionosphere, and since the latter rotates together with the earth, so too do the 
field lines higher up and together with them the whole ion-electron plasma. The 
ions are thus compelled to maintain a constant angular velocity (rather than an- 
gular momentum) and the resulting centrifugal force further increases the scale 
height. The net result is that at a height of one earth radius protons become more 
numerous than neutral hydrogen, their concentration decreasing very slowly with 
height. We thus find that the inner radiation belt is populated primarily by neutral 
and ionized thermal hydrogen and the outer belt primarily by thermal protons and 


electrons. 


The presence of this ionized hydrogen plasma, or protonosphere as Johnson?25 
has called it, surrounding the earth to large radial distances was discoveree3® and 
investigated by means of “whistlers’’. In general the ionospheric layers are a reflect- 
ing medium for all frequencies below about 5 Mc/s. However, at very low (audio) 
frequencies the ionosphere is transparent to the axtraordinary component of electro- 
magnetic waves which travel in the direction of magnetic field lines. Further, the 
field lines have a facussing actior on the wave which will follow the magnetic lines 
from one hemisphere to the other and then back again. Now electromagnetic energy 
in the audio region is propagated by lightning discharges, and by measuring the 
time delay of each frequency component to make the round trip along field lines 
from hemisphere to hemisphere, the electron density at the top of the mafenetic 
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path may be estimated. In this fashion Storey36 found unexpectedly high electron 
densities of several hundred per cm3 at distances of several earth radii, and since 
then the total distribution has been estimated by this method37. 


It may be mentioned that the magnetic field line through Haifa reaches a height 
of 2000 km above the magnetic equator, and it is hoped that whistler measure- 
ments being conducted by the lonospheric Physics group of the Technion will 
yield information of thermal ion concentration at the base of the inner radiation 
belt and its correlation with electron densities at the F layer peak (as obtained by 
ionosonde techniques). 


The total ion content of the ionosphere has been calculated by measuring the 
Faraday effect (the rotation of the plane of polarisation of an electromagnetic wave 
in a magnetoionic medium) for radio frequency signals transmitted from satellites3®. 
These measurements have verified that 2 to 3 times as many electrons may be found 
above the F layer peak as below. 


CONCLUSION 


The nature and generation of the Van Allen radiations, the thermal whistler medium 
or protonosphere, and the telluric neutral hydrogen corona have been described. 
A great deal concerning the role of the solar wind and the extended solar corona 
is not clear as yet, and their clarification awaits the more extensive measurements of 
the future earth satellites. 
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BUCKLING OF THIN CIRCULAR CONICAL SHELLS SUBJECTED 
TO AXISYMMETRICAL TEMPERATURE DISTRIBUTIONS* 


JOSEF SINGER 


Department of Aeronautical Engineering, Technion-Israel Institute 
of Technology, Haifa 


ABSTRACT 


A method developed previously for the analysis of the instability of conical shells 
under external pressure! is now applied to the problem of buckling due to axisym- 
metrical temperature distributions. 

The method is based on a solution of the Donnell type stability equations, derived 
by Seide? ia the presence of slightly relaxed boundary conditions. The solution 
satisfies the usual simple support conditions regarding the radial deflections w, 
rigorously, but implies elastic restraints as far as the axial and circumferential 
displacements are concerned. However, the constraint on the circumferential, v, 
displacements is of such a nature that it practically represents the rigid fixation 
required by the usual simple supports; and the much weaker constraint on the axial, 
u, displacements (in the direction of the generators) combines with the very small u 
displacements themselves to a condition approximating freedom from axial restraint. 
The effect of the u restraint on the instability of cylindrical and conical shells under 
external pressure was previously investigated and the critical pressures found to 
differ by 1-2 percent for typical shells*°. 

The stability equation in the radial direction is then rederived in a modified form 
to facilitate solution by the Galerkin method. 

This solution may be applied to any type of loading for which the prebuckling 
stresses are functions of x only and the hoop stresses are the prime cause of buckling. 
The thermal stresses due to an axisymmetrical temperature distribution are such a 
case. 

The thermal stresses may be computed from Meissner type equations, since the 
equilibrium equations upon which the Seide stability equations are based, reduce 
to Meissner equations for rotational symmetry. A solution for the Meissner equations 
in the presence of axisymmetrical temperature distribution given by Huth®, but ina 
slightly modified form, is used here for the calculation of the thermal stresses. 

For the thermal buckling analysis the thermal stresses are then expressed so that 
their magnitude depends on a temperature parameter 7T,, which also determines the 


* This paper is based on research supported by the United States Air Force through the Air 
Force Office of Scientific Research, ARDC, under Contract No. AF 61 (052)—339. 
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magnitude of the temperature distribution, and their functional form is approximated 


by power series in x. 

A method for the calculation of the critical temperature parameter T, (which 
indicates the onset of buckling), similar to the above mentioned method for the 
buckling under external pressure, is then developed. The method is also suitable 
for analysis of a shell subjected simultaneously to external pressure and axisymmet- 
rical temperature distributions. Typical cases are analysed. 
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ON THE PHYSICAL RIEMANN-CHRISTOFFEL TENSOR 
IN ORTHOGONAL COORDINATES 


Zv1 KARNI 
Division of Mechanics, Technion-Israel Institute of Technology, Haifa 


ABSTRACT 


Expression for the physical Riemann-Christoffel tensor is derived and its geometrical 
meaning stated. 


1. The derivation of the physical Riemann-Christoffel tensor, also known as the 
curvature tensor, is based on the concept of the curvilinear derivative of a physical 
tensor and that of the physical Christoffel symbols. The definition and derivation of 
these concepts form the subject of a separate paper* of which only the main results 
will here be stated. 

A three-dimensional curvilinear network is formed by the families of curves x; = 
= const. (i = 1,2,3) where the x; denote the curvilinear coordinates of a point in the 
space. The network can also be represented by means of the position vector r(x,,X2,X3) 
or r (yi, ¥2, ¥3) where the y, stand for the cartesian coordinates. The positive unit 
tangent vectors to the curves at a point of the network, denoted t,, form a local 
cartesian system of coordinates. A physical tensor is defined with respect to the 
loca] system, hence it is always a cartesian tensor. 

The two-fold metric system g,; (which is not a physical tensor) is found from 


_ [Ya PY (1.1) 
ij 6x: Ox 


where it should be understood that repeated Greek indices only imply summation 
while repeated Latin indices are not summed. The necessary and sufficient condition 
for a curvilinear network to be orthogonal is 


Jij =9 (i#j) (1.2) 


The discussion in this paper is confined to orthogonal coordinates only, hence all 
indices are written as subscripts. 


~ * By Braun, Karni and Reiner, to be published shortly. 
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The ‘‘curvilinear derivative’? at a point P of the network of a vector field A is 
identified with the vector difference of the vector field while proceeding along a curve 
of the network and per unit arc length ds;. This derivative is denoted by a comma, thus 


OA <3 
eile (1.3) 
0s; 
where 
ds, = ./ 94 dx; (1.4) 


The physical Christoffel symbols G;;,, are defined as the projections of the three 
principal curvatures 1/p; of the coordinate curves on the planes of the local system t; 
at P. Thus G;,, means the projection of the principal curvature 1/p; on the plane 
formed by the unit tangents t; and t,, namely the plane t; x t,. From the definition 
it follows that 

Gizx = =Gi; (1.5) 


Furthermore, for three mutually different indices 
Gijxn=0 (i AJ#K Fi) (1.6) 


because the principal curvature 1/p; has no projection on its tangent t;=t; x t. 
In view of (1.5) and (1.6), the matrix of the non-vanishing symbols takes the form 


[Gil] = — Gu, = 


The G;;, can be calculated from the components of the metric system g;; by means 
of the relation 


ee eee uence 
2 Gi ey du 


1 nes ey a) (1.8) 


| Gn G32 | (1.7) | 


It should be noted that the physical Christoffel symbols form a tensor of the third | 


rank with the property of skew-symmetry — Eq. (1.5). 
The i-th component of the curvilinear derivative of A—Eq. (1.3) —is now expressed 
in terms of the symbols G;;, as 


aA, (1.9) | 


(A,,); = 4,3 = ae (ALG ys 


a) 


=-% 
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The curvilinear derivative of an absolute physical tensor of rank n is given by 


Biered: 
A; i, +i ft A a Ay et, Opn oe Ay 


6 Oya isto 
na? ads jai,, ( ) 


i;i2° oe 


2. The physical Riemann-Christoffel tensor R;j is obtained on performing the 
difference in the second, mixed derivative along a curve of the vector A;, namely 


Aita— Aijue = Aa Rien (2.1) 
In view of Eqs. (1.8), (1.9), 1.10) and the identity, obtained from Eq. (1.8), 


Giix + G ix: + Gyij — 0 
it is found that 


Rija = a = - + GygiGipi — Gigi Gegj; + Goi; Gorr (2.2) 
Again, in view of Eqs. (1.5), (1.6), (1.8) it is readily established that 
Rijut = — Rijn =e R jixi co Ryij (2.3) 


The symmetrical and skew-symmetrical properties of the physical tensor are thus 
similar to those of the corresponding analytical tensor, denoted here as Grip 
which is obtained on performing the similar difference in the second covariant deri- 
vatives of a contravariant vector A‘. The principal distinction between the two 
tensors is in the dimensions of the two. All the compononts of the physical tensor 
possess the same dimension, that of the square of curvature, whereas the components 
of the analytical tensor need not be of the same dimension in the various curvilinear 
coordinates. The relation between the two tensors in orthogonal coordinates is 


Ria = = 6/99 Gadu Rijn (2.4) 


The six non vanishing components of the Riemann-Christoffel tensor in a three 
dimensional space ,as results from Eq. (2.3), form a symmetric matrix 


| Ryat2 Riz2z-Riasi | | 
|| Riga ll 1.23 =} - Ry323 Ra331 | (2.5) 
; : R3131 | 
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It is also known that 
R; 
oe eaten en 
911922 


where K stands for the Gaussian curvature; hence from Eq. (2.4) 
Riaz = —K (2.6) 


We therefore conclude that the three ‘‘normal’’ (diagonal) components in the matrix 
(2.5), namely of the type R;j;;, coincide precisely with the Gaussian curvatures of — 
the three surfaces formed by each pair of the coordinate curves of the network. 


3. If we consider the vector t, the components of which ft; are the unit tangents | 
of the network at P, as an example of a vector field, then the single non-vanishing | 
component of the Riemann-Christoffel tensor in a two-dimensional space (surface) | 
is found from Eq. (2.1) and (2.3) to be 


ty 12—ti,21 = t2Ry212 = Ri212 = t2,21 — F212 = 1 Roi21 = Raia (3.1) 


Now, t;,12 means the projection on the local axis t, of the second vector difference | 
of the tangent t while traversing along the arc of unit length ds, fiirst, then along 
the unit arc ds,, namely 


aaa (3.2) 


: 
OS OS 4 


Alternatively, if the tangent vectors at P,P’,P’,P” are labelled t’ t’, t’, t” respectively 
(Figure 1b), then t, ,. is obtained from the vector difference of these vectors taken in 
the order 


warzaly 


Similarly, t; 5, is obtained from the vector difference taken in the order 


(t”—t’)—(t’ —t) 
or CCE | 
thay == (er oe (3.3) | 
OS; OSaJ 


These projections, on either t; or tj, are in general not equal and the Riemann- 
Christoffel tensor measures the difference in value of these projections. In the case 
of the sphere r(0,y,a) for example (Figure 1a) 
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) Figure 1 
| a. Curvilinear network on a sphere of unit arcs to ds; 
b. Tangent vectors with respect to the local systems 
the matrix of the physical Christoffel symbols reads 
| 0 | 
IGulls2=| ete? 0 | (3.4) 


| a 
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making use of Eqs. (1.1), (1.3), (1.9), (1.10) and (3.4), we find for Eqs. (3.2), (3.3) 


ctg 0\2. 
12 47 i) 


a 


10 /ctgé 
—ty 21 ='o,21 = lee 


Consequently, 


lo ;sctgé ctg 0\2 1 
Rya12 = Raia = ae =n) ele SKS 


aé0\ a a a2 
in view of Eq. (2.6). 


What has been said about the geometrical interpretation of a single component 
of the physical Riemann-Christoffel tensor applies equally well to the other compo- 
nents, all being of the same physical dimension. Not so with the analytical tensor 
the meaning of which can only be sought in conjunction with other quantities which 
make the entire expression an invariant. We know that if a vector A“ of unit length 
undergoes a parallel displacement from a point P along a closed curve in a con- 
tinuum, then the angle formed by the vector on returning to P can be expressed in 


terms of R;,, and the surface element bounded by the curve. When, therefore, the 


direct physical or geometrical interpretation of a tensor is of an interest, the 
physical tensor has this advantage over the corresponding analytical one. 


HYDRODYNAMICAL APPLICATIONS OF A THEOREM ON SPHERICAL 
MEANS 


Ert JABOTINSKY 
Division of Mathmatics, Technion-Israel Institute of Technology, Haifa 


ABSTRACT 


A Theorem on vectorial spherical means is applied to the angular momentum theorem 
for a spherical volume completely immersed in a laminar flow. The resulting formula yields 
directly classical results for perfect fluids and shows that rotretroty is an indicator of 
magnitude of the moments of shear and thus possibly also of the conditions for the 
inception of turbulence. 


Consider a portion of the medium ina volume t inside a closed surface o. Let dt be 
an element of the volume t. Denote by p(x,y,z,t) v(x,y,z,t) and a(x,y,z,t) respec- 
tively the density, the velocity and the acceleration of the matter in dt. We note that 


Let fod be the outside volume force acting on the matter in dt, and let Fdo be the 
outside surface force acting through do. Furthermore denote byr the position vector 
of points in t or on a, that is the position vector of volume elements dt or of surface 
elements do relatively to a fixed origin O. 

Then the angular momentum theorem, applied to the matter within a, yields: 


a, (r x F)do + ph (rx f)pdc = “fff x v)pdt. 


Differentiating under the integral sign and using the equation of continuity an 


the fact that = F and collecting the volume integrals, we find: 
[£@x Edo = Sf fe x (a feat. (1) 


We shall need here only the special case in which (1) is applied to the matter inside 
a spherical surface © and the origin O of the position vectors r is the fixed point 
which, at the moment t, coincides with the center of £. Here only the tangential 
components of the surface forces F contribute to the moments (rx F)do. We shall 
denote by S the spherical volume inside £. 
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THEOREM I. 
Consider a continuous medium of density p(x,y,z,t). Inside that medium con- 
sider a spherical surface & of radius R, enclosing the spherical volume S. Let r 
be the position vector of the surface element d& of X or of the volume elements dS 
of S relative to the center of X. Let TdX be the shear applied to the matter in S by 
the outside medium through d& and let pi(x,y,z,t)dS be the external volume force 
applied to the matter in the volume element dS. Then, if a is the acceleration of 

the matter in dS: 
4 


[J @x Td = F2R* rot p(a— Ff), (2) 


where rot p(a — f) is the weighted average of rot p(a —f) in S given by: 


a Je 5 rot p(a — f)dS = 7g rot pla —f). (3) 


Proor: The angular momentum theorem can here be written 
Sfo@x Dade = ff f,fex@- fleas, 
and it remains to transform the volume integral. We have: 


[ff (x @—Aleds = f ap [Jy, ex (@- Mody, 


where XZ, is the spherical surface of center O and radius f. Putting r=fn where n 
is the unit vector normal to dZ, and denoting by S, the interior of Zs, we have, by 
a Classical result: 


Jf x (a — f) odd, = BJS {n x (a — f) odd, = BS, rot p(a — f) dS, 


so that: 


Nhe [x x (a — f)]pdS = [pap J Toe — f)dS, = ‘CEE —- r* rot pa —f)dS. 


We thus have: 


[fox Dak = fff <>" rot p(a—fyas, 4 


and the main theorem follows because: 


Rie te, ae 
(GH) 5 aS = GR’. 
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APPLICATION TO SOLIDS 


In a solid rota is constant in space. Indeed if @ is the angular velocity of the solid 


This is due to the fact (which is easily checked analytically) that in a solid: 


and: 
rotv = 2a. 


Therefore, for a homogenous solid, our Theorem I becomes: 
THEOREM II. 


(Homogenous solid medium). If the medium is solid and homogenous of cons- 
tant density p then: 


Sffy@ x Tat = ar °p . (6) 


eye 2 ‘ ‘ ; : d@ ‘ ee 
This yields interesting conclusions on the maximum ap solid can sustain without 
splitting. 


APPLICATION TO PERFECT FLUIDS 
In a perfect fluid T = 0 and Theorem I yields easily: 
THEOREM III. 


(Perfect fluids). In a non-viscous medium where all shears are zero: 
rot p(a—f) = 0. (7) 


Three special cases are of interest: 

A) Incompressible perfect fluid. 

Here p is constant. If f is the constant field of gravity, rot f= 0 and Theorem IV 
reduces to: 
THEOREM IV. 

In an incompressible perfect fluid under gravity: 


rota =0 (8) 


B) Two dimensional horizontal flow of a perfect fluid. 
Vector a is now horizontal while f is vertical. Writing that the vertical component 
of rot p(a — f) is zero, we find, in this case: 
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THEOREM V. 


In a two dimensional horizontal flow of a perfect fluid: 


rot pa = 0. (9) 


C) Two-dimensional horizontal irrotational flow of a perfect fluid: 
In this case rotv = 0, which implies rota=O (this is easily checked analytically). 
Remembering that 
rot pa = prota—axXgradp and applying formula (9) we find: 
aX grad p = 0 (10) 
which yields: 


THEOREM VI. 


In a two-dimensional irrotational flow of a perfect fluid in a_ horizontal 
plane the acceleration is parallel to the gradient of the density. 


APPLICATION TO VISCOUS FLUIDS 


We limit ourselves to the incompressible case with constant viscosity uw. The Navie- 
Stokes equation here is: 
' p(a—f)= — gradp + wAv 


But rot grad p = 0 and Av = grad div v — rot rot vy = — rot rot v for div v = 0, 
so that Theorem I now becomes: 


THEOREM VII. 


In the flow of an incompressible fluid of constant viscosity: 


R2— r2 4n 5. = +e 
Lf x T)d= = — uff. ; —rotrotrotydS = — oe Lrotrotroty 

This result shows that the value of | rot rot roty| is a measure of the moments of 
the shears in the considered flow and may possibly be anindicator for the condition | 
of the inception of turbulence, 


THE MECHANICAL PROPERTIES AND STRUCTURE OF SOME 
HIGH TENSILE ALLOY STEELS 


S. BECKER AND M. G. BADER 
Israel Atomic Energy Commission and Battersea College cf Technology, London 


ABSTRACT 


Some outstanding results have been obtained while the mechanical and tempering properties 
of some silicon - molybdenum — vanadium — copper high tensile alloy steels have been in 
vestigated, in the water quenched and tempered condition. Metallography has been limited 
to optical microscopy, in spite of this limitation some relationship between mechanical 
properties and microstructure could be found. 


1. PREFACE 


The work described in this report was carried out at Battersea College of Technology, 
London, and forms part of a much larger investigation in progress. 

During the period of the present werk the first author has been in receipt of a 
Travelling Scholarship provided by the British Council, for which he would like to 
express his appreciation. 

The authors would also like to thank all those in the Metallurgy Department, 
Battersea College of Technolcgy, for their advice and encouragement, especially 
to Mr. L. W. Derry and Dr. B. L. Daniel. 


2. MATERIAL 


The steels were manufactured from Swedish iron, and appropriate master alloys 
in a vacuum induction furnace. Ingots of approximately 56 Ib. weight were produced 
and these were hot rolled to a bar of approximately 5/8” diameter except in the case 
of steel “D’” which was rolled to approximately 24” x 5/8” so that test pieces 
could be extracted in the transverse direction. Test piece blanks about 23" long 
were heat treated and test pieces were then machined. 
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The composition and treatment of the steels are shown in Tables I and II and it 
should be noted that Steels “A” and “‘D” are of similar composition and treatment 
but in “A” the test pieces are parallel to the rolling direction and in “D” they are 
transverse. 


On account of the high hardness of many samples of both, the gauge length and 
the threaded portions of the test pieces were ground. 


Owing to the short supply of material there were in some cases insufficient test 
pieces for a full range of heat treatments. When this was so, small test specimens 
were heat treated to obtain the hardness figures. 


TABLE I 
Analysis of steels (%— WT) 


Steel C Si Cu Mo Vv Mn S P Ni Cr Al 


C 0.45 P22 0.90 1.86 0.13 0.35 0.019 0.060 0.02 0.04 0.015 
Q 0.38 1.56 2335 0.98 0.35 0.47 0.009 0.005 0.01 0.02 0.04 

) eee OER I/ 1.50 1.62 2.08 0.36 0.52 0.008 0.026 0.01 0.04 0.04 

OF 70°32 1.50 IRA 1.49 0.34 0.43 0.009 0.018 0.01 0.03 0.04 
E 
F 
A 


0.29 (hye 2.06 0.94 0.28 0.50 0.027 0.025 0.02 0.04 — 
0.28 1.06 1.45 2.50 0.34 0.25 0.023 0.034 0.02 0.04 — 
0.18 1.50 d73, 0.80 0.30 0.29 0.033 0.016 0.06 O16 — 
DS OLS 1.50 ATS 0.80 0.30 0.29 0.033 0.016 0.06 0.16 — 


Ae = Aluminium: 


TABLE II 
Heat treatment 
All specimens quenched from 970°C and tempered to the following tenperatures: 


Spec. Treatment Treatment Treatment Treatment Treatment Treatment Treatment Treatment 


No. Cc Q IE, O E F A D 
1 — A.Q. _ A.Q. A.Q. A.Q. A.Q. A.Q. 
Zz 87°C 100°C 100°C 100°C 105°C 110°C 100°C 100°C 
3 215°C 245°C 235°C 235°C 205°C 220°C 200°C 200°C 
+ 395°C 39556 400°C 400°C 390°C 405°C 390°C 390°C 
5 505°C $35°C 515°C $25°E 515°C 530°C 505°C 505°C 
6 600°C 600°C 600°C 600°C 600°C 600°C 600°C 600°C 
i — — 650°C 650°C 650°C 650°C 650°C 650°C 
8 — —- — 700°C 700°C 700°C 700°C 700°C 
A.Q. = as water quenched. 


The 0.1 per cent proof stress, ultimate tensile strength, elongation reduction of area, 


Young’s modulus and hardness have been measured in duplicate experiments and 
plotted graphically. 


More than 50 microsections have been examined and 15 photomicrographs are 
presented. 
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3. EXPERIMENTAL PROCEDURE 


(a) Mechanical testing 


One hundred and eleven test pieces were tested on a Denison 50 ton tensile testing 
machine, using a 1” Lamb optical extensometer. Full stress/strain graphs were 
plotted and elongation and reduction of area of test piece was measured. Typical 
stress strain curves are reproduced in Figures 1-8.* 

Hardness tests were conducted on all specimens on a Vickers diamond pyramid 
hardness testing machine with 30 kg load. Microhardness tests were conducted on 
some of the specimens using a Bell Telephone testing machine. 


(b) Metallographic examination 


Some fifty microsections cutin the longitudinal and transverse direction from 
the tensile test pieces have been examined after being polished and etched with2% 
Nital. Twenty-three photomicrographs illustrating typical structures are shown. 
Only 10 photomicrographs are shown herewith. 


4. MECHANICAL TEST RESULTS 


The results of the tensile tests are shown in Figures 9-17.** 
The following quantities were determined: 
(a) 0.1 per cent proof stress (Tons per sq. inch) 
(b) Ultimate tensile strength (Tons per sq. inch) 
(c) Elongation (% on 4, / Area) 
(d) Reduction of area (%) 
(e) Hardness (V.P.N. HD/30) 
(f) Young’s Modulus (Ib/sq. inch) 
Steel “C” (Figures 1 and 9) 
This steel has the highest carbon content (0.45%C) of those tested but the other 
alloying elements are relatively low, particularly copper and vanadium. 


A tensile strength of about 150 tons per sq. inch. together with elongation of 8% 
was obtained in the as quenched condition. The tensile strength drops qui e rapidly 
on tempering to 100 tons per sq. inch at 600°C. (This may be compared with 0.45 
per cent plain carbon steel as quenched approx. 130 t.s.i. dropping to 45 t.s.i. on 
tempering at 600°C). 

The hardness follows the same trend as the tensile strength and agrees with the 
usually quoted relationship of V.P.N. x 0.2 = U.T.S. Proof stress values however 
show an increase with tempering with a slight drop as the higher tempering tempe- 


* Only Figures, 1,6,7 are included in this paper (pp. 70-71). 
** Only Figures 9,10,13,14,15,17 are shown herewith (pp. 71-74). 
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ratures are achieved. With steel “C’” as with most of the others investigated the 
elongation stays nearly constant on tempering up to 500°C while the reduction of 
area increases initially and then decreases. 


Steel O 

This steel is of lower carbon content than “‘C” and its alloying elements are some- 
what higher. With this steel the rate of decreasing tensile strength and hardness is 
less than in “‘C”. The 0.1 per cent proof stress increases on tempering up to 600°C, 
the ductility increases on tempering and the ratio of proof stress over ultimate tensile 
strength approaches unity. 


Steel “L” 

There is doubt whether all the specimens inthis batch responded to heat treatment. 
In particular there is a 'arge difference between test pieces 2 and 2A*. As it was not 
possible to obtain more material, small specimens cut from the broken test pieces 
were heat treated and the hardness figures obtained were very different from those 
obtained from the original heat treatment. 


This being so, no further evaluation can be made on any of the other properties 
of the material or its structure. 


Steel “O” 


The tensile strength shows only slight drop on tempering with a tendency to secon- 
dary hardening around 600°C. The proof stress rises slightly and then more dras- 
tically at 600°C, above 600°C the strength drops rapidly. 


Steel“ E” 


This steel is very similar to steel “‘O”’ which is discussed above. 


Steel ‘‘F’’ (Figures 6 and 14) 

Generally very similar to Steels “O” and “E” but the secondary hardening is 
more marked. On tempering at 600°C its tensile strength is equal to the as quen 
chéd condition. This steel has the highest amount of carbide forming elements (Mo,V). 


Steel “A” (Figures 7 and 15) 

This steel has the lowest carbon content (0.18 %C), of all the steels tested. 

The general pattern of the curves is similar to steel ‘“‘F”. The tensile strength and 
hardness after tempering at 600°C are even higher than in the as quenched condition, 


also the 0.1 per cent proof stress shows the greatest increase between the as quenched 
and 600°C tempering temperatures of the steels under investigation. 


** Two specimens representative of each composition and heat treatment have been tested. 
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Young’s Modulus 


This follows the well known!3 trend i.e. decrease of modulus for qunenched steels 
with increase of carbon content. From Figure 17 it is seen that Young’s modulus of 
steels “C” and “Q” is lower than of steels “L”, “O”, “BR”, “F’, “A” which also 
show a tendency of increase of modulus on tempering. The modulus also increases 
quite appreciably on secondary hardening. 


5. METALLOGRAPHIC EXAMINATION 


Low power examination of longitudinal sections of the broken tensile test pieces 
revealed considerable banding in most of the steels. This can be seen in plates 1-4. 
There was also a considerable quantity of non-metallic inclusions often associated 
with the banding but not uniform throughout the bars. 

Isolated large areas of a light etching constituent occurred in some steels, particu- 
larly in those of lower carbon contents as shown in plate 1. The exact nature of this 
constituent is not understood but it is possibly ferrite. The micro-hardness is 200-220 
D.P.N. 

It is possible that complete homogenisation was not obtained during the melting 
process and that these are regions of almost unalloyed iron. 

In conclusion, it can be seen from the low power observation that the steels are 
far from homogenous and also relatively dirty for high grade alloy steels. This may 
be explained by the small experimental melts and small amount of reduction in hot 
working to produce the test bars. The need to conserve material has probably resulted 
in more of the ingot being utilized than in normal steelworks practice. 

Examination at high power was carried out on several sections from each steel. 
In the as quenched condition the higher carbon steels “C’’, ““Q”’, ““L” and “O” were 
predominantly martensitic and the lower carbon steels “E”, “F” and “A” showed 
considerable amounts of bainite. : 

Fine particles of microconstituent were visible in most of the as quenched steels 
and in some cases delineated the previous anastenite grain boundaries. On tempering 
they increased in size and number, and in some cases the grain boundary effect was 
most marked. (Plates 6, 16, 20, 23). Precipitation also occurred in a random manner. 
The precipitated particles were very small and it was not possible to establish whether 
they were carbides or possibly the ¢ (Cu Rich) Phase. 

In most cases the background structure changed only slightly on tempering, and 
the sorbitic structure was only apparent on tempering at 650-700°C. 


6. DISCUSSION 


The breakdown of martensitein quenched steel during tempering has been studied 
by a number of workers, in carbon steels using mainly X-ray techniques, and dilato- 
metric methods. In recent years with the development of the electron microscope it 
has been possible to follow the constitutional changes during different stages of 
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tempering. Using these means it has been proved that quenched hardened low alloy 
steels have the structure of tetragonal martensite, they may also contain residual face 
centred cubic austenite and “‘auto-tempered martensite’ 5>11+13, 

The change in physical properties during tempering results from the combined 
effect of the decomposition of the two solid solutions, martensite and anstenite which 
is observed in three stages: 

(i) Precipitation of ¢ - iron carbide (below 150°C) 

(ii) Transformation of retained anstenite into bainite (150°C-280°C) 

(iii) Precipitation of cementite (above 200°C) 

(iv) In some alloy steels an alloy carbide may form on tempering usually in the 
400°C-700°C range. 

The metallography in this work has been limited to optical microscopy which does 
not allow a detailed observation of the tempering process to be made. In spite of 
this limitation an attempt will be made to find a relationship between the mechanical 
results with: 

(a) Composition 

(b) Heat treatment 

(c) Microstructure. 

Generally different quenching methods and alloy additions may alter the nucleation 
conditions of carbides or of solid solution and precipitation hardening. In designing 
for alloy steels to be used in the quenched and tempered condition it is usual to assume 
that if efficient quenching is achieved the as quenched strength of the steel depends 
predominantly upon the carbon content. The alloying elements are effective mainly 
in controlling martensitic hardenability and tempering characteristics so that in 
alloy steels a considerable upward displacement of the temperature ranges of tempe- 
ring processes can occur. 


It has been shown that when plain carbon steels are tempered from martensitic 
condition within the temperature range of 400°C-700°C their hardness decreases 
logarithmically with tempering time?. If the steels contain non-carbide forming 
alloying elements, this general type of softening still prevails although the rates of 
softening are decreased (e.g. silicon 0.5% delays martensite breakdown into ferrite and 
cementite due to solid solution hardening 10, similar considerations apply to copper). 
With the addition of strong carbide forming elements such as molybdenum and 
vanadium the softening behaviour or tempering is greatly modified by the formation 
of a dispersion of alloy carbides. Between 200°C and 400°C i.e. within the temperature 
range of the precipitation and growth of cementite, alloy elements retard softening 
caused by tempering (see references 5» 13) when present in sufficient amounts these 
alloy elements are about equally effective in retarding softening or even increasing 
hardness. If the amount is not sufficient to combine with all the carbon the retarding 
effect is not so obvious. 
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Thus there is a rapid decrease in hardness in a plain carbon steel between 250° and 
450°C, but in the steels under investigation this drop has been delayed until tem- 
peratures between 600°C and 700°C are reached. 


The addition of copper apart from retarding tempering of martensite confers some 
solution hardening by dissolving in ferrite and when present in sufficient quantity 
precipitation hardening occurs by precipitation of the ¢ phase (solution of Fe in 
F.C.C. copper), this latter effect has by far the greatest influence on strength. The 
maximum hardening occurs in the range 400-500°C in copper bearing plain carbon 
steels 4, but the temperature and extent of hardening might well be influenced by 
other alloying elements present. 


Development of present steels 


In the development of the steels investigated there have been two main aims: 


(a) Reduction on volume change on quenching, associated with the austenite to 
martensite reaction, and thus reduce internal stress and possibility of cracking on 
quenching. 

(b) To retain high strength even after tempering at high temperatures. 

To achieve the above aims silicon and molybdenum have been chosen as the prin- 
cipal alloying elements as they reduce the volume change during transformation 
10, 12, 14 to a greater extent than more conventional alloying elements. In addition 
both elements have a marked effect in reducing martensitic breakdown on tempering. 

Resistance to softening in the medium temperature range is achieved by the pre- 
cipitation hardening due to the copper, and at higher temperatures secondary harde- 
ning occurs due to the formation of molybdenum and vanadium carbides. 

The full interpretation of this somewhat complicated state of affairs is beyond the 
scope of the present work but, nevertheless, certain trends are indicated by study of 
the data obtained: 

(a) Quenched strength is determined by carbon content. Steel “Q” bas approxi- 
mately the same composition as “E” except for the carbon content, steel “Q” showg 


a higher quenched hardness than steel “E”’. 

(b) The copper to carbon and the silicon to carbon ratios determine initial decrease 
in strength on tempering 9 10. 

This trend is seen mainly when comparing steel “Q” and “A” which differ mainly 
in the carbon content other alloying elements being similar in both steels. [t is sug- 
gested that the higher the ratio of S; + C,/C the less wiil be the decrease in hardness 
in the 200°C-450°C tempering range. 

(c) The ratio of molybdenum and vanadium to carbon determines secondary 
hardening. Steel “F” having twice the ratio of steel “C” will show more pronounced 
secondary hardening. This may indicate that there is a certain ratio 0° Mo + V versus 
carbon for secondary hardening. This was also observed by some investigators in 


other steels 2) 5. 11, 
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(d) The limit of proportionality is very low in the as quenched steels as can be 
seen in Figures 1, 6, 7 and in the P.S./U.T.S. ratio. On tempering there is an increase 
in P.S./U.T.S. which approximates to 0.90 at 600°C. 


This is typical of precipitation hardening alloys. It is noteworthy that an increase 
in limit of proportionality occurs even on tempreing below 200°C corresponding 
to the first stage of martensite breakdown (eé carbide formation), or to transformation 
of retained anstenite although there is no metallographic evidence of this. 


7. CONCLUSIONS 


Remarkable mechanical properties are obtained in the steels investigated. Tensile 
strength of 90-110 tons per sq. inch together with elongations of the order of 10-15% 
are obtained with little or no reduction on tempering up to 600°C. 


These steels are still in the experimental stage and have yet to be evaluated in ser- 
vice but would seem to offer improved properties when compared with more conven- 
tional steels, particularly for operation at elevated temperatures. 


In their present form water quenching is necessary for hardening and while it is 
more usual to prefer oil quenching steels in industry the objection to water quenching 
is largely overcome, because these steels have a low volume change on transformation 


There is insufficient evidence from the present work to suggest an optimum compo- 
sition for this class of steel, and it is probable that variations for different applica- 
tions will be required. However, although higher strength is obtained in steel ““C” 
by nature of its higher carbon content the ductility is reduced, which suggests that 
the lower carbon versions might be more suitable for general use. 


It would no doubt be possible to produce an oil: hardening version of these steels 
by increasing the hardenability by addition of suitable elements. 
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PRELIMINARY DEVELOPMENT OF AN ANNULAR JET INJECTOR 


A. KOGAN AND M. VICTOR 
Department of Aeronautical Engineering, Technion-Israel Institute of Technology, Haifa 


ABSTRACT 


The paper describes research carried out at the Department of Aeronautical Engineering 
on an injector designed for operating the Department’s transonic wind tunnel. Experiments 
were carried out in the low and high-subsonic ranges with a view to determining the in- 
fluence of the injector components on its performance, using mixing tubes of variable 
length and two different diffusers. 

It was established that optimum mixing tube length ranges from Lm/Dm= 2.5 to 4 
and that a 4° diffuser permits improvement of performance. It was also found that the 
performance level of the injector in question is lower than those of similar systems described 
in references. This may be attributed to the negative influence of the transition piece between 
the rectangular and circular sections, separation of the secondary flow before the injection 
area and the geometry of the injection slot. 

Suggestions for improvement of performance are given. 


1. INTRODUCTION 


The experimental transonic wind tunnel (Figure 1) operates on the induction prin- 

ciple, consisting in the injection of high-pressure air (inducing air, or primary flow) 

into a channel in the direction of flow required for the induced air (secondary flow) 
As a result of the mixing process, momentum is imparted to the induced air and the 

desired flow is obtained. Injection may be central or peripheral (Figure 2); in most 
cases the latter method yields better results. 


In wind tunnels, injection is effected downstream of the test section in order to 
- preclude interference with the secondary flow. The smallest cross section of the 
inducing air nozzle (injection slot) is normally located in the injection area. Although 
a sonic velocity is obtainable in the slot when the stagnation pressure of the primary 
flow (blowing pressure) is about double the static pressure in the injection section, 
it is usually advantageous to use a higher blowing pressure. While this would not 
produce a heigher velocity in the slot, it increases the density of the primary flow 
and hence the momentum imparted to the induced air. The same consideration 
holds if the minimum cross section is located before the injection area. 


Injector performance may be evaluated by measuring the ratios p,/p,; of the 
blowing pressure and secondary-flow stagnation pressure required for obtaining a 
given flow at Mach number M, in the test section and m,/m, of the respective rates 
of mass flow. The performance level rises with increasing mass ratio and, conversely 
with decreasing pressure ratio for a given mass ratio and Mach number. 
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Figure 1 
Induction-type transonic wind tunnel. 
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The suitability of the induction system for the operation of wind tunnels was 
studied in a number of laboratories!"8. In the low subsonic range, the blowing 
pressure required is low and the mass ratio very high, but the induction method is 
impracticable in view of the large storage capacity required for the inducing air at 
low pressures. In the supersonic range, above M > 1.5 in the test section, the blow- 
ing pressure is very high and the mass ratio decreases rapidly, so that the induction 
method is impracticable in this case as well. In the transonic range, however, the 
method was found practicable. The blowing pressure required for obtaining transonic 
velocities is not too high and the mass ratio range is reasonable. 


2. DESCRIPTION OF INJECTOR AND HIGH-PRESSURE AIR SUPPLY SYSTEM 


The injector in question is of circular cross section and injection section diameter 
of 110 mm. The advantage of the circular cross section lies in the simplicity cf its 
design3>4, The test chamber, i.e. the minimum effective cross section area (see §5), 
is rectangular, dimensions 69 x 89 mm. A transition piece from rectangular to 
circular cross section is interposed between the chamber and the injector. The in- 
jector, designed for open-circuit experiments, is shown in Figure 3 with a conical 
diffuser attached at its exit. 


The supply of inducing air to the pressure chamber is effected from tanks through 
valve V and a 1” pipe. A throat p* is fitted in the pipe for regulating the inducing 
air mass. After the throat, the pipe branches out providing two 1” inlets for the 
high-pressure air (Figure 4). 

2.1 Construction 


The air inlet A (design following Ref. 9) is made of sheet steel, and has a 
convergence ratio of 40:1. The transition piece B is made of aluminium casting, 
with a 70 x 90 mm end attached to the air inlet. These are also the dimensions of 
the test chamber to be interposed between the air inlet and the transition piece. 
These dimensions ensure condensation—free flow in the transonic range}%1!!. This 
part comprises a gradual transition from rectangular to circular cross section of 
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Figure 3 
Injector 


100 mm diameter. In order to preclude choking of the secondary flow in the trans- 
ition piece? or near the injection slot>>7, a 1.24 area-ratio was adopted for the circular 
and rectangular cross sectional areas. Since this part is in effect a diffuser, its total 
opening angle does not exceed 6.5°. 


Cis a fixed steel tube of 100 mm internal diameter and E, screwed into D, a similar 
tube of 110 mm internal diameter, this variation of diameter being an additional 
precaution against undesirable choking of the secondary flow3-7. The gap between 
parts C and E comprises the variable injection slot (Detail A, Figure 1). The slot 
is designed for continuous adjustment by rotating E with respect to C with the aid 
of handles M. The shape and dimensions of the slot provide for a 0-3 mm width 
range against a 25 mm horizontal displacement of E, ensuring sonic velocity in the 
slot and a constant injection angle of 8° between slot and injector axes. 


Since in the case of peripheral injection perfect slot symmetry is essential, special 
attention had to be paid to precision in manufacture of the slot with respect to the 
concentricity of parts C, D and E. The total angle of conical transition from the 100 
mm diameter of part C to the 110 mm of part E was 14°. 
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The primary flow enters the annular pressure chamber D, designed for a blowin g 
pressure of about 40 atm, through pipes P, and then part E through the injection 
slot. The secondary flow follows a path through the air inlet A, transition piece B 
and tube C, and meets the primary flow ina mixing tube consisting of parts E and F 
(Figure 4). This mixing tuube has a constant diameter of 110 mm and its length is 
adjusted by varying the attachment F, providing a range of length-to-diameter 
ratios between 4 and 10. The F pieces are fitted with conical diffusers G with the 
required angle and area expansion ratio. 


3. MEASUREMENTS AND INSTRUMENTATION 
Evaluation of injector performance requires the following data: 

Mach number in the test section, blowing pressure, primary and secondary rates 
of discharge and cross sectional area of the slot, as well as the static pressure distri- 
bution along the injector and stagnation pressure distribution at the mixing tube 
(or diffuser) exit. 


3.1 Velocity measurement in the test chamber 
The Mach number in the test chamber was determined by measuring the mean 
static pressure at the four walls of the chamber (pressure taps i, Figure 4). Assuming 
that the stagnation pressure in the chamber equals the atmospheric pressure 


teen! 
Z Dae wy? i 
sia era, 1] ©) 
Static and atmospheric pressures were measured by means of a Collesman manometer. 


3.2 Blowing pressure measurement 
Since the ratio of the pressure chamber and slot areas is 70 for a maximum slot 
opening and increases with decreasing slot area, the static pressure measured at the 
chamber wall closely approximates the stagnation pressure. This pressure is read 
by means of a Bourdon gauge. 


3.3. Primary discharge measurement 

The primary rate of discharge can be measured by two methods: (a) by determining 
the stagnation pressure, throat area and stagnation temperature in the high-pressure 
air pipes (Figure 4), assuming that the flow in the throat is sonic; or (b) by deter- 
mining the same data in the pressure chamber on the same assumption for the slot 
In both cases the rate of discharge is given by Eq. (2), assuming further that flow 
from the stagnation point is isentropic and the discharge coefficient (k* or k; res- 
pectively) is 1. 
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Figure 4 
Injector for transonic induction-type wind tunnel. Location of pressure measurement taps and inducing air supply. 
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In the first case the throat area (A>) is known and constant. The pressure measured 
before the throat (p,,) is not the true stagnation pressure since at that point the 
velocity is not zero. The throat was proportioned so as to maintain sonic flow even 
at the minimum slot opening. The pressure measured at pressure tap,, (Figure 4) 
is 0.994 of the stagnation pressure and readings should be corrected accordingly 
At pressure tap, the static pressure of the inducing air is measured after the throat 
The ratio of this static pressure and the stagnation pressure before the throat gives 
a rough indication of sonic flow in the throat. The mass ratio can be obtained ac- 
cordingly with the aid of Eq. (2), for all cases where the PplPop Tatio is below the 
critical level (0.5283). 


The stagnation temperature before the throat cannot be measured directly. It 
may be assumed, however, that it lies between those of the atmosphere (i.e. of the 
air in the tanks) and of the pressure chamber, which can be measured. The difference 
between these two is fairly small (10°C), so that it can be assumed that the stagnation 
temperature of the inducing air before the throat equals the atmospheric. 


In the second case, the inducing air mass can be calculated with the aid of Eq. (2) 
if the blowing pressure is about double the atmospheric, ensuring sonic flow in the 
slot. The slot area (A,,) can be determined by measuring its width at minimum 
opening with the aid of standard-thickness shims and counting the number of half- 
turns of part E (Figure 4). The minimum area is zDs,. The slot width increment 
due to the displacement of part Eis As = n/2 x 1/16 sin 8°, where nis the number of 
half-turns (the number of threads per inch being 16 and the injection angle 8°, 
Figures 1 and 4). Hence, for a given n, the slot width is: 


S=S + As (3) 


and the slot area: 
Ajg = xDs (4) 


3.4 Secondary discharge measurement 


The induced air mass is given by 


A; M; 
of aioe? vee bain 


(1 Z zie 


p,; and T,,; being atmospheric pressure and temperature resp., and the discharge 
coefficient k; = 1. 
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3.5 Determination of mass ratio 


The mass ratio is given by 
ivan 
(—) aa Pot As é = Bere M; (6) 
A 


Pie * +1 
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calculation being based on the stagnation pressure in the throat. 


Alternatively, the mass ratio calculation can be based on the blowing pressure y he 


pr. 
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3.6. Measurement of static pressure distribution along injector and stagnation pressure 
distribution at exit 


For measuring the static pressures, taps were provided along the transition piece, 
mixing tube and diffuser. The stagnation pressures at the mixing tube or diffuser 
exit were measured with the aid of pressure rakes. All measurements were carried 
out with the aid of an alcohol multimanometer or mercury manometers, according 
to the range of pressure. 


4. PRELIMINARY TESTS 


The first test series was carried out with mixing tube No. 2 (length-to-diameter 
ratio 6) fitted with a conical diffuser with a total opening angle of 8° and area ex- 


pansion ratio 8, designed in accordance with experimental results given in Ref. 
2525) 5 and 6. 


4.1 Location of test chamber 

The test chamber is defined as the zone of lowest static pressure of the secondary 
flow. In order to locate this section, a series of pressure taps was provided at the end 
of the air inlet and along the transition piece (Figure 4). The injector was operated, 
for two slot widths at three different velocities of secondary flow. This velocity was 
determined for each width by adjusting the blowing pressure by means of valve V. 
Results (Figure 5) show that (a) the section of the lowest static pressure is about 
60 mm from the transition piece inlet for all velocities of secondary flow; (b) flow 
in the transition piece is steady; (c) slot width has no influence on pressure distribu- 
tion in the transition piece. This last indicates that the injection zone is sufficiently 
removed from the test chamber so that at subsonic velocities flow in the latter 
unaffected by the injection process. 
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Figure 5 


Pressure distribution along the transition piece as function of injection slot area and test section 
Mach number. 


With the test chamber located as described above, four pressure taps were provided 
on its walls for measuring the mean static pressure. Separate measurements showed 
a maximum discrepancy of about 0.1” mercury between readings, i.e. within the 
limits of experimental error (accuracy + 0.05” mercury). 

The Mach number in the chamber is determined by the average pressure p, in 
this section. 


5. EVALUATION OF INJECTOR PERFORMANCE 


5.1 Procedure 


For a given slot width determined by the number of half-turns of part E, the inducing 
air valve is opened so as to produce a static pressure corresponding to a desired 
value of Mach number M,, see Eq. (1). Readings were taken for n = 0, 0.5, 1, 2, 
3, 4, 5, 7, 10, 15, 20, 25, and for a Mach number range from 0.4 to 1 at intervals 
of 0-1: 

5.2 Determination of ‘“‘aerodynamic” slot area 


As explained in §3.3, one of the principal factors governing injector performance 
is the slot section area. Measurement of the initial width sy by means of standard 
shims entails a relative error ranging from 20% forn =0 to 5% for n = 15. Since 
the slot area is a parameter of the peformance curve!7, the latter cannot be plotted 
unless this error is minimised. One means of achieving this is “aerodynamic” mea- 
surement of the section area. It can be assumed with fairly close accuracy that flow 
in the inducing air supply pipe is one-dimensional with normal shock wave, involv- 
ing no other losses between the throat ,.« and the slot ;. Heace, with sonic flow in 
both throat and slot, we have 


Disp” oa Poj Aj (8) 
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For each given n, corresponding stagnation pressures before the throat and the 
slot were obtained for different Mach numbers in the test chamber. Since the throat 
area is known, the slot area can be calculated for each case with the aid of Eq. (9) 


A, [sq.in.] = 0.1028 222 (9) 
oj 


in which A¥ = 0.1028 sq.in. is the throat area. 


The mean “aerodynamic” area for each n, together with the corresponding mean 
deviation, is given in Table I. 


TABLE I. 
Mean “aerodynamic” area for mixing-tube of NALS aii ratio 6, 8°-diffuser of area expansion 
ratio 8. 
n Aj 61 x 103 6:/Aj 
Half-turns **Aerodynamic” Mean yA 
of part E mean area deviation 
*0 0.170 — — 
*.0.5 0.199 — — 
1 0.228 P49 1a3 
2 0.289 6.45 2.10 
3 0.360 3.95 1.09 
4 0.424 4.22 1.00 
5 0.490 3.50 0.71 
i 0.619 1.90 0.31 
10 0.806 9.25 215 
15 1.118 9.00 0.80 
20 1.426 Sit 2.20 
25 1.730 — — 


* Calculated by extrapolation from Figure 6. 


The area, calculated by means of direct geometric measurement as explained in 
§3.3, is plotted against n in Figure 6, giving the result of the preliminary series. 
The diagram shows that the ‘‘aerodynamic’”’ value (A,) is always higher than the 
one obtained geometrically (A,,), the maximum deviation not exceeding pe AS 
hence the error entailed in direct measurement is systematic and the mean “aero- 
dynamic” area can be assumed as the actual slot area. 


Since for small openings (n < 1) the flow in the throat was not sonic, as expected 
(see §2.1 and §3.3), A, was not calculated with the aid of Eq. (9), but determined 
by extrapolation of the curve in Figure 6. The final “aerodynamic” area was de- 
termined as the mean of the complete series of experiments. It is presented in Table II. 
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Figure 6 
Injection slot area. Results based on geometric and aerodynamic calculations. Measured initial slot 
width Sy = 0.098 in; 8° diffuser. 


TABLE II. 


Mean “aerodynamic” area, based on complete series 


n 
Half-turns 0 0.5 1 2 3 4 5 Wy 10 15 20 DS) 
of part E 
A ,, sq. in. 0.172 *0.201 0.225 0.295 0.364 0.430 0.495 0,621 0.808 1.121 1.435 1.730 
AJA 0.018 0.021 0.024 0.031 0.038 0.045 0.052 0.065 0.085 0.118 0.150 0.182 


* Calculated by extrapolation. 


5.3. Determination of pressure ratio 


In all tests the stagnation pressure of the induced air was equal to the atmospheric. 
Results given in Figure 7, show the pressure ratio for the original injector (mixing 
tube length-to-diameter ratio 6; 8° diffuser). It is about 100% higher than that of 
the best similar system reported in the references?»3 for small slot openings, and 
about 50% higher for the larger openings (Figure 8). For simplicity, the diagram is 
confined to the Mach number extremes of 0.4 and 1.0, but the same applies for the 
entire range. 
5.4 Determination of mass ratio 

The mass ratio was calculated with the aid of Eq. (6) and (7). Results (Figure 9) 
show complete agreement between the two methods of calculation for A,/A, = 0.04 
and M, > 0.5. For lower Mach numbers and smaller openings, the discharge co- 
efficient k,, not allowed for in Eq. (7), deviates consid. rablly from 15. It was also 
found that for 4,/A, < 0.024 flow in the throat is not choked, so that results based 
on Eq. (6) are no longer accurate for this range. 
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Figure 7 


Blowing pressure as function of slot area and test section Mach number. Injector with 8° diffuser 
of area expansion ratio 8 and mixintg tube of length-to-diameter ratio 6.0. 


Figure 10 shows part of the results given in Figure 9 compared with corresponding 
results given in Ref. 2. For small openings, the mass ratio obtained in the present 
experiments is 50% lower than the corresponding ratio in Ref. 2, and for the larger 
openings—25% lower. Since the overall limit of accuracy does not exceed 6%, 
this discrepancy cannot be attributed to limitations of the measurement technique 
but to increasing losses in the system in question. 


6. INFLUENCE OF INJECTOR COMPONENTS 


6.1 Transition piece 
A divergent transition piece causes an average reduction of about 10% in perfor- 
mance level, in both the subsonic and supersonic ranges2:3, With the transition 
piece used in the present experiments higher losses should not be expected, since 
in the range in question there is no separation of flow (Figure 5). 
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Blo wing pressure as function of slot area and test section Mach number for different injectors. 


6.2 Divergent section between part C and injection secticn 
Table III shows the static pressure in section a before the slot (Figure 4) for different 
Mach numbers. The effective diameter of this section was calculated with the aid 
of one-dimensional Flow equations. The values obtained are 5% to 10% lower 
than the geometric diameter of 110 mm, indicating separation of the flow before 
this section. Considering that the total opening angle of this section is 14° (see 


TABLE II 
Effective diameter of secondary flow in injection section 


Pa/Patm pen pair 8, Aa Da 
M; Mean of Ma = mm2 mm 
complete series A* Ai Ai/A* 
0.4 0.948 0.277 2.19 1.38 8480 104.0 
0.7 0.871 0.448 1.45 132. 8100 101.5 


1.0 0.815 0.548 1.26 1.26 7750 99.4 


92 A. KOGAN AND M. VICTOR Bull. Res. Counc. of Israel 


tf 


RL 
Se 
Wi 


q | at fos, 


—#- MASS RATIO BASED ON Ay (equetion7 ) 
—e— MASS RATIO BASED ON Ap (equation 6) 


2 i 


| | 


Co) 0.02 0.04 0.06 0.08 0,40 0.42 
INJECTION SLOT OPENING RATIO &y/Ai 


Figure 9 


Mass ratio as function of slot area and test section Mach number. Injector with 8° diffus er of area 
expansion ratio 8 and mixing tube of length-to-diameter ratio 6.0 


§2.1), it can be assumed that for a smaller angle, performance will be improved 
considerably. 


6.3 Injection slot 
6.3.1 Injection angle 
Experimental results3, show that this angle has almost no effect on performance, if 
smaller than 10°. This includes the present series, in which the angle was 8°. 


6.3.2 Shape of slot 

The slot, as shown in Figure 1, comprises a long stretch of constant cross section. 
For small openings the length-to-width ratio is 40, so that a pressure drop can be 
expected. This drop, calculated with the aid of one-dimensional flow equations 
allowing for friction, ranges from 1% to 6% according to opening. Minimising 


this loss necessitates a modified shape as shown in Figure 11, with a reduced ef- 
fective length-to-width ratio. 
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Figure 10 
Mass ratio as function of slot area and test section Mach number for different injectors. 
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Figure 11 
Cross section of proposed injection slot 
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Figure 12 


Peripheral distribution of actual injection slot width, as determined by geometric measurement. 
Injection slot opening ratio 0.018 (n = 0) 


6.3.3 Non-symmetry 


The stagnation pressure distribution at the exits of the various mixing tubes (Figure 
13) and diffuser (Figure 14) show that the flow is non-symmetrical with respect to 
the injector axis. This non-symmetry is especially pronounced for small slot openings, 
irrespective of blowing pressure. 

The width of the slot was measured at different points along its circumference by 
means of standard shims, showing a deviation from symmetry of about 0.2 % (Figure 
12). The effect of this non-symmetry on injector performance was studied in Refs. 
2 and 3. The loss in blowing pressure due to this 0.2°% deviation should not exceed 
6% for an opening ratio of about 0.04. 


6.4 Mixing tube 


The longer the mixing tube, the more complete the mixing process, and the higher 
the momentum imparted to the induced air. Above a certain optimum length, how- 
ever, the momentum increment is smaller than the loss due to increased friction. 


In order to determine the optimal length of mixing tube, a series of experiments 
was performed with the mixing tubes listed in Table IV. 


In all tests the static pressure distribution was measured along the tube, with the 
diffuser attached and without it. Where no diffuser was used, the pressure was found 
equal to the atmospheric near the point x,,/D,, = 4.5; it rose further along an ad- 
ditional stretch of the tube and eventually dropped again to atmospheric level (Fig- 
ure 15). The corresponding mass ratios are given in Figure 20. 


With the diffuser attached, the maximum pressure point was shifted backwards 
to about x,,/D,, = 2.5 for an 8° diffuser (Figure 17) and to Xml Dm = 4 for a 4° 
diffuser (Figure 19). The corresponding blowing pressures and mass ratios are 
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Stagnation pressure distribution at diffuser exit along vertical diameter for Mi = 1.0. 4° diffuser of 
area expansion ratio 8. 


TABLE IV 
Mixing tube data 


Serial Length Length to 
No. Lin (Cm) diameter ratio 
Lm/Dm 

1 44 4 

2 66 6 

3 88 8 

4 110 10 

Di 63 Siy// 

33 80 3 

yY 99 9.0 


given in Figures 21 and 22, respectively. It can be concluded that, for the injector 


in question, the optimum mixing-tube length-to-diameter ratio ranges from 2.5 
to 4. 


The problem of optimum mixing-tube length was also studied in Refs. 5 to 6. 
According to Ref. 5, the optimum ratio for an annular injection system without a 
diffuser is about 5, while according to Ref. 6, the optimum ratio for a central system 
with a diffuser is 6.25. The fact that the former requires a shorter tube is due to the 
larger initial contact surface between inducing and induced air. It is assumed that 
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Figure 15 
Static pressure distributions along mixing tubes. Pg; = 20 psig. Injector without diffuser. 


the irregularity of the secondary flow before the injection section, due to its separa- 
tion in the transition piece between C and E (see §6.2) is a contributing factor to 
the acceleration of the mixing process in the experiments described here. 


6.5 Diffuser 

The static pressure distribution along the 8° diffuser with an area expansion ratio 
of 8 (Figure 16) indicates separation of the flow in all tests. This is also indicated by 
the stagnation pressure distribution at the diffuser exit. The flow separation zone 
(Figure 18) coincides with the section segment in which the air velocity at the mixing 
tube exit is lower (Figure 13). 

In Figure 17, the static pressure distributions for both diffusers attached to mix- 
ing tube No. 2, are compared with the pressure distribution calculated by one-di- 
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Figure 16A 
Static pressure distribution along mixing tube No. 1 and 8° diffuser of area expansion ratio 8. 


mensional incompressible flow theory. This comparison permits location of the 
flow separation zone for the diffuser with the 8° opening angle. No appreciable 
separation was detected in any of the 4° diffuser tests. 


The 4° diffuser permitted an average 12% improvement in pressure and mass 
ratios (Figures 8, 10, 23 and 24, Table V). 


7. SUMMARY AND CONCLUSIONS 


The conclusions drawn from the tests described above are as follows: An increase 
in mixing-tube length above L,,/D,, = 2.5 or 4 respectively leads, in all cases, to a 
reduction in performance level. Elimination of the flow separation in the diffuser 
by reducing its opening angle results in improved performance. Pressure distribution 
along the mixing-tubes and diffusers is unaffected by slot width. 


Vol. 9C, 1961 A. KOGAN AND M. VICTOR 99 


M.=6.4 


—— Aj/Are 0.034 


+ Aj/ai = 0448 


Xn [Dyn | 
0 o4 O2 O3 OF OF O68 Q7 O8 08 40 


Ka /be 
xine | 
DIFFUSER at rah 


TUBE 


Figure 16B 
Static pressure distribution along mixing tube No. 2 and 8° diffuser of area expansion ratio 8 


The best performance level in the present tests was found to be 60% below 
those reported in Refs. 2 and 3 for low slot opening ratios. Even with the higher 
opening ratios, a difference of about 20% was still found. The main cause is 
believed to be separation of the secondary flow before the injection section. It is 
expected that a reduced opening angle in part C will lead to considerable 


improvement. 


Additional loss factors are the presence of the transition piece between the rect- 
angular and circular sections, imperfect symmetry and high primary flow friction 
in the slot. In future tests within the framework of this study, the system will be 
modified in the light of these conclusions. It is also proposed to study the effect of 
the diffuser inlet shape and area expansion ratio. 
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Figure 16C 
Static pressure distribution along mixing tube No. 3 and 8° diffuser of area expansion ratio 8. 


TABLE V 


Injector performance with 4° and 8° diffusers compared with Ref. 2. Mixing tube length-to-diameter 
ratios 6, Mach number in test section 1.0. 


Poj Poj Poj Poi Poi Poj 
AjlAi psia psia psia ons es Pa, 73 (Poi) ner 
Ref. 2 4° diffuser 8° diffuser 4° diffuser 8° diffuser 
0.02 100 160 195 1 1.60 1.95 
0.07 40 57 61 1 1.43 1.60 


0.12 35 42 46 1 1.20 1.34 
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Figure 16D 
Static pressure distribution along mixing tube No. 4 and 8° diffuser of area expansion ratio 8. 
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— gas constant 

— injection slot width 
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— specific heat ratio 

— absolute error (Table 1) 
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Static pressure distribution along diffusers compared with theoretical. Injector with Mixing tube 
No. 2 and Aj/Ai = 0.118. 


PLOW DIRECTION 


Figure 18 
Zone of separated flow in the 8° diffuser as determined from stagnation pressure measurements 
at diffuser exit. 
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Figure 19C 


Static pressure distribution along mixing tube No. 3’ and 4° 
diffuser of area expansion ratio 8. 
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Figure 19D 


Static pressure distribution along mixing tube No. 4’ and 4° 
diffuser expansion ratio 8. 
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Figure 20 


Mass ratio as function of 
mixing tube length and 
injection slot area. Injector 
without diffuser. 
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Blowing pressure and mass ratio as functions of mixing tube 
length. Injector with 8° diffuser of area expansion ratio 8. 
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Figure 22 
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Blowing pressure and mass ratio as functions of mixing tube length. Injector with 4° diffuser of 


area expansion ratio 8. 
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Blowing pressure as function of slot area and test 

section Mach number. Injector with 4° diffuser of 

area expansion ratio 8 and mixing tube of length-to- 
diameter ratio 6.0. 
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Figure 24 


Mass ratio as function of slot area and test section Mach 
number. Injector with 4° diffuser of area expansion 
ratio 8 and mixing tube of length-to-diameter ratio 6. 


AERODYNAMIC RESEARCH FACILITIES AT THE DEPARTMENT OF 
AERONAUTICAL ENGINEERING, TECHNION-ISRAEL INSTITUTE 
OF TECHNOLOGY 


JOsEF ROM (RABINOWICZ) 
Technion Israel Institute of Technology, Haifa 


ABSTRACT 


The research facilities developed by the staff of the Aeronautical Engineering Department, 
Technion-Israel Institute of Technology, are described. The aerodynamic facilities include 
the 60 x 90cm closed return subsonic tunnel, and the 60 x 60 cm open circuit low 
speed tunnel. Two blow down supersonic tunnels 10 x 10 cm and 30 x 30cm fed bya 
high pressure air supply of a 1000 ft. at 200 atm. An induction transonic tunnel model 
7 x 9cmanda 7 meters long 7.5 x 7.5 cm shock tube. The report describes briefly 
these facilities andtheir associated instrumentation. 


I. INTRODUCTION 


Experimental aerodynamic research is based on the various facilities for simulating 
the flow conditions, as well as instrumentation for the measurement and recording 
of the experimental results. 

The aerodynamic laboratory at the Technion was established in 1955. The cons- 
truction of the high pressure air supply was completed in 1957 and the 10 x 10 cm 
supersonic tunnel was completed in the beginning of 1958. The 60 x 90 cm closed 
return subsonic tunnel was completed at the end of 1958. During 1959 — 1960 the 
open circuit low speed tunnel was added. Design and construction of the 30 x 30 cm 
blow down supersonic tunnel was undertaken. A transonic tunnel model driven by 
an air injector was designed and built, and a shock tube of 7.5 x 7.5 cm square 
section 7 meters long is under construction and will be completed early in 1961. 

Concurrently, design and development of various instrumentation systems was 
undertaken. A number of manometers, multimanometers and bi-fluid micromano- 
meters were installed. A three component side balance for the subsonic tunnel was 
designed and and built. Under development now are various sting balances for all 
aerodynamic force and moment components. The recording system include, strain 
gauge bridges, pen-recorders and oscillographs, as well as, optical schlieren system 
for the supersonic tunnels. 

The present report comes to summarize the development and to describe briefly 
these facilities which are used in the research and teaching activities of the Aero- 


nautical-Engineering Department. 
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II. SUBSONIC WIND — TUNNELS 


A. The subsonic 60 x 90 cm closed return wind tunnel 


The subsonic tunnel is a conventional closed return tunnel. The design methods and 
overall characteristics of this type of tunnel are described in Ref. 1. An overall view 
and major dimensions are shown in Figure 1. The tunnel is constructed of a wooden 
frame covered by a sandwich structure of plywood and mazonite, except, for the 
propellor section which is housed in a steel cylinder. The tunnel is powered by a 20 
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REMARMs DIMENSIONS IW PARANTHESES. DESCRIBE DeetH ,oF SECTION: 


Fiy 1 + The Go Xo0™ Subsonic Wind Tunnel 


H.P. a.c. motor provided with a mechanical speed reduction system. The power is 
transmitted by belts to the propellor. The maximum propellor rotation speed is 1000 
RPM. The propellor is mounted on a central body and is follwed by straightner 
vanes. The air circuit includes then the return duct, corners, settling chamber, con- 
traction section, test section, diffuser and returns via two more corners to the propeller 
section. The corners are equipped with circular vanes in order to reduce the losses. 
A turbulence reducing screen of high solidity is installed in the settling chamber. 
The contraction cone has a 6 to | contraction ratio. 


The test section has a 60 x 90 cm cross section and is 100 cm. long. Both top, 
bottom and side walls are removable. The side walls are made of transperent perspex 
enabling a full view of the test section. At the end of the test section there is a 1/4” 
wide atmospheric vent. 


The main diffuser has 7° total opening angle. A safety screen for model catching 
is provided at the end of the main diffuser. 
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All tunnel cross sections are rectangular in shape except the round propellor 
section. Appropriate transition section from rectangular to round are provided 
both fore and aft of this section. 


The tunnel top speed is 50 m/sec (110 mph) and the turbulence level is below 1 %- 
Velocity distribution are uniform over the test section area within 1/2 %. 


B. Open circuit 60 x 60 cm tunne! 


A simple open circuit tunnel has been constructed using a ventilation blower powered 
by a 7.5 HP. d.c. motor which is controlled by a motor generator set. The blower is 
fitted with a central cone and the flow is directed through a transition section to a 
square channel of 60 x 60 cm cross section. This channel is constructed of 50 cm 
long sections. The overall length can be adjusted according to the testing demands. 
Various interchangeable contraction sections have been built. These are then connec- 
ted to the end of the square channel so that a number of test sections of various cross 
sectional area and air velocity can be used. A short diffuser is also added to the 
test sections in order to increase the tunnel’s efficiency and air spead. Figure 2 shows a 
typical arrangement of this tunnel. The versatility of this tunnel makes it particularly 
suitable for student experiments and other programs. 
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III]. TRANSONIC AND SUPERSONIC FACILITIES 


A. High pressure air system 
An air storage of 1000 cu.ft. at 200 atm. is used as a prime-mover for a number of 
blowdown high speed wind tunnels, and other facilities. 

Three Worthington 4 stage reciprocating compressors each one powered by a 
100 H.P. a.c. electrical motor are used in parallel to charge up the high pressure 
containers. The compressed air is cooled and the oil and moister are separated by a 
flow reversal stage and then by a circumferential vortex separator. The cooled and 
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dried air is then pipped into the steel storage bottles. Two 2 1/4” high pressure pipes 
bring the air to a 3” pneumatically controlled ““Conoflow” valve. The outlet of the 
control velve can be directed by various valves to the blow-down facilities. Safety 
diaphragms are fitted in the circuit to ensure relief in case of a blockage at some point 
in the low pressure side of the system. A diagram of the system is shown in Figure 3. 


B. Supersonic blow down wind tunnels 


A 10 x 10 cm blow down supersonic wind tunnel is in operation. This tunnel has 
the following structural sections:- (a) settling chamber, (b) nozzle and test section 
(c) diffuser (d) exhaust pipe. The settling chamber is a cylindrical tank 18” in diameter 
and 40” long. The air is admitted through a 90° angle directed backwards so that 
good mixing is enhanced. The forward flange is provided with an entrance cone 
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Figure 4b 
10 cm Supersonic Blowdown Wind Tunnel 
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which has a transition to a rectangular section of 4” x 7’. The nozzle section is a 
rectangular box 4” wide 7” high and 26” long, in which the nozzle blocks are attached 
to the top and bottom walls and so provide a test section of 10 x 10cm. The boundary 
layer correction is obtained by increasing the test section size in the vertical plane 
only. A 6” window is opened to the test section for optical observations. Nozzle 
blocks for Mach Numbers of 1.5; 2.25 and 3.5 are available. The diffuser is housed 
in a similar box. Wooden liners are used for defining the diffuser sections. This 
diffuser has a fixed second throat with an opening slightly larger than the best the- 
oretical area. A model holder with adjustable angle of attack range is incorporated 
in the diffuser section. This mechanism is based on the parallelogram method in 
which the model holding bar follows the distortions of an external parallelogram 
about an axis which is centered in the test section. 


From.the diffuser section the air is led through a 12” exhaust pipe to the atmos- 
phere. An overall view of the tunnel is shown in Figure 4. 


A second branch of the air supply is directed to the 30 x 10 cm tunnel. It is similarin 


overall layout to the loxiocm one with all linear dimensions multiplied by approxi- 
mately 3. 
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C. The7 x 9 induction transonic wind tunnel model 


A small induction tunnel which is used as a model for a transonic facility is now 
under development. An air injector supplied by the laboratory high pressure air 
supply is used for inducing a large mass of air through a 7 x 9 cm test section. 
The tunnel is designed for operation at Mach numbers between 0.7 to 1.4. 

The tunnel has a slotted wall test section, a circumferential injector followed by 
a lenght of mixing tube and a conical diffuser. The tunnel is shown in Figure 5 
and is described in Ref. 2. 


D. Shock-tube 


A shock tube is now under construction in the laboratory. The shock tube is driven 
by either compressed air or hydrogen at pressures up to 3000 psi. The compression 
chamber is a thick walled 3” steel pipe 1.5 m long. A metal diaphragm separates 
the high pressure from the low pressure sections. A square section 7.5 x 7.5 cm 
and 7 m long can be evacuated to less than 1 mm Hg abs. The low pressure sec- 
tion is divided into a number of sections of which the two sections comprising the 
last 3m, are provided with instrumentation ports and observation windows. The 
low pressure chamber ends in a dump tank. High capacity vacuum pump is con- 
nected to the dump tank. A diagram of the tube is shown in Figure 6. A descrip- 
tion of the shock tube and its instrumentation is given in Ref. 3. The | m section 
at the end of the low pressure chamber can be replaced by a hypersonic nozzle. 
This facility will be used for research in high energy, high Mach number 


flows. 
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IV. INSTRUMENTATION 


The various laboratory facilities impose a wide range of requirements on the in- 
strumentation systems to be used. In general the systems can be catalogued accor- 
ding to the measured end recorded quantity. In most aerodynamic research the | 
physical parameters which are of interest are the pressure, temperature, density, 
forces and moments and in supersonic flows also the shock waves position and 
shape. A block diagram describing the various building blocks of an aerodynamic 
instrumentation system is shown in Figure 7. 


A. Pressure and temperature measurements 


Pressure is measured on model surfaces or in the flow by special pressure rakes. 
Some pressure rakes used in the 60 x 90 cm subsonic tunnel are shown in Figure 8. 
The pressure from the orifice at the model surface or in the rake is transferred through 
special tubing to manometers or other pressure gauges. Number of multimano- 
meters, with fixed and varying inclination of the tubes, using colored alcohol are | 
available. Special high sensitivity bi-fluid micro-manometers using kerosine-al- 

cohol were developed. For higher pressure mercury manometers are also available- 


The manometers and the bourdon type pressure gauges are read and recorded 
manually or photographed. Records are then tabulated and evaluated. 
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Figure 8a 
Pressure Rake for the 60 « 90 cm Subsonic Wind Tunnel 


Figure 8b 
Total Pressure Rake for the 60 ~ 90 cm Subsonic Wind Tunnel 


The short running times in the blowdown supersonic tunnels and the 1ecording 
of transient pressures is obtained by the use of strain gauges pressure transducers. 


A number of Statham unbonded strain gauge transducers for pressure ranges bet- 
ween 0.5 to 500 psi are available. The transducer output is recorded cn a pen re- 
corder or an oscillograph. Ovtput can be also displayed on an oscilloscope when 
required. 

Temperatures are obtained by various thermocouple combinations Cu-Cons- 
tantan or Iron constantan and output is recorded on a self balanced potentiometer 
or an oscillograph. 
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B. Force and moments measurements 


The restriction of the model sizes for the various facilites require small force ba- 
lances. Development of side balances and sting balances was undertaken. A three 
component side balance using strain gauge instrumented flextures is used in the 
60 <x 90 cm subsonic wind tunnel. This balance can also be adapted to the super- 
sonic wind tunnel. The balance is shown in Figure 9. 


The internal sting balances have bakelite strain gauges which are bonded to the 
balance structure and respond to all componenents except drag. The drag element 
is measured by the deformation of a flexture arrangement which is sensed by an 
unbounded strain gauge. These balances are described in Ref. 4 and are shown 
in Figure 10a and 10b. 


The strain gauge balance output is recorded on a pen recorder, an oscillograph 
or an oscilloscope. 


The sting balances are calibrated on a special test stand (Figure 11). All six force 
and moment components can be loaded simultaneously in order to evaluate all 
interacticns. 


Figure 9 


Three Component Side Balance 


Figure 10b. Five Component Drag Balance 
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AN APPROXIMATE METHOD OF SOLID PROPELLANT INTERNAL 
STAR GRAIN CALCULATION 


A. YARON AND R, CoRETT 
Scientific Department, Ministry of Defense 


ABSTRACT 


The usual internal burning grain configurations are calculated with the sole con- 
sideration of the burning perimeter variation. The proposed method concerns itself 
mainly with the overall dimensions of a grain and permits the general outlay of an 
internal star as well as the main proportions of the charge. The method of calculation 
takes into consideration the total mass of the charge, the ballistic properties of the 
propellant and permits calculation of the approximate shape of the internal star, 
which may then be modified by the usual methods according to required variations 
of the burning perimeter. The proposed method provides an easy means of charge 
calculation giving a possibility of quick appreciation of the motor proportions. 
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ON EXTRAGALACTIC COSMIC RAYS 


H. KASHA 


Department of Physics, Technion-Israel Institute of Technology, Haifa 


ABSTRACT 


It has been established by experiments 1:2 on the extensive air showers of cosmic 
radiation that there is no evident change in the character of the energy spectrum of 
primary cosmic rays in the energy range E > 1016 ev, and that particles of at least 
1019 ev are present in the primary cosmic-ray flux 13. These facts impose stringent 
conditions on the theories of origin of cosmic rays, which in general assume that 
the particles are stored and mixed (or even accelerated) by the galaciic magnetic 
clouds. The efficiency of this process is impaired at ~ 10!7 ev, which should lead 
to a deficiency and a pronounced anisotropy of the high-energy taii of the radiation. 
The absence of such phenomena is difficult to reconcile with a purely galaciic theory 
of origin of cosmic rays. In his “supernova” theory, Ginzburg+ completely neglects 
any extragalactic contribution. This may be inconsistent with this very successful 
theory itself; many galaxies, in addition to our own, possess cosmic-ray sources 
(supernovae) and some special extragalactic objects like the spherical galaxy M87 
are sources of bremsstrahlung radiation (and thus presumably cosmic rays) on an 
incomparably larger scale than the galactic supernovaes-6, The contribution of 
extragalactic cosmic-ray sources depends greatly on the “transmission coefficient” 
of the galactic boundary’. The value of this coefficient can be, at present, only guessed 
at; however, the experimental facts on air showers at very high energies seem to 
contradict clearly some estimates, in particular that of Ginzburg, who put it at ~ 1%. 
The transmission coefficient is certainly energy-dependent (increasing with increasing 
energy), and thus it is not very convincing to exclude the possibility of an extraga- 
lactic origin of at least a part of the cosmic-ray particles initiating the greatest exten- 
sive air showers. 


REFERENCES 


CLARK, EARL, KRAUSHAAR, LINSLEY, ROSSI AND SHERB, 1957, Nature, 180, 353, 406. 
ALLAN, BEAMISH, BRYANT, KASHA AND WILLS, 1960, Proc. Phys. Soc., 76, 1. 
ZATSEPIN, G. T., 1960, Proc. Moscow Cosmic Ray Conf., vol. I. 

GINzBurG, V. L., 1958, Progress in Cosmic Ray Physics, 4, 337. 

BAADE, W., 1956, Astrophys. J., 123, 550. 

SHKLOovskil, I. S., 1955, Astron. J. Soviet Union, 32, 215. 

Si1TE, K., 1956, Bull. Res. Counc. of Israel, 6A, 61. 


NDAAPYN 


120 Bull. Res. Counc. of Israel, Vol. 9C, 1961 


STRAIN GAUGE BALANCE FOR WIND TUNNELS 


A. KOGAN, J. SCHWARTZ AND A. SEGINER 
Department of Aeronautical Engineering, Technion Israel Institute of Technology, Haifa 


ABSTRACT 


Two electrical resistance strain gauge balances were developed for wind tunnel 
model testing: one—for the measurement of lift, side force, pitching moment, yawing 
moment and rolline moment, and the other—for the measurement of lift, drag, side 
force, pitching moment and yawing moment. 


An effort has been made to reduce the interactions between the different force 
and moment components to a minimum. 


A special unit has been developed for the drag measurement. It consists of a 
system of unbonded strain gauges. 
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